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Abstract 

Neutrinos are the only particles in the Standard Model of particle physics that have only 
been observed with left handed chirality to date. If right handed neutrinos exist, they could 
be responsible for several phenomena that have no explanation within the Standard Model, 
including neutrino oscillations, the baryon asymmetry of the universe, dark matter and dark 
radiation. After a pedagogical introduction, we review recent progress in the phenomenology 
of right handed neutrinos. We in particular discuss the mass ranges suggested by hints for 
neutrino oscillation anomalies and dark radiation (eV), sterile neutrino dark matter scenarios 
(keV) and experimentally testable theories of baryogenesis (GeV to TeV) . We summarize con- 
straints from theoretical considerations, laboratory experiments, astrophysics and cosmology 
for each of these. 
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This review is intended to give a comprehensive overview of what we know (more precisely: 
what I know) about the phenomenology of right handed neutrinos. 

Faced with the difficulty of writing a text that is readable within reasonable time, and is at 
the same time is as precise as possible, I provide extra details, explanations and comments for 
the interested reader in an extensive set of footnotes. For a quick read you may essentially ignore 
all of them. 

Furthermore, each section has been written to be as self-contained as possible, so that they 
can also be read individually as "mini reviews" . The basic notations used throughout the review 
are introduced in sections 11.11 and [21 A very quick summary is given in appendix El 

I will maintain a list of corrections online sX^http://res-puhlica.eu/RightHandedNeutrmos.html\ 
If you find any errors or typos I would be very grateful if you could inform me; you can find my 
most up-to-date contact information on the above website. 
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1 Introduction 



The Standard Model of particle physics (SM) and theory of General Relativity (GR) form the basic 
pillars of modern physics. Together they can describe almost all phenomena observed in natur^ 
in terms of a small number of underlying principles - general covariance, gauge invariance and 



quantum mechanics - and a handful of numberqj All elementary particles we have observed 
to date can be understood as fundamental excitations of a few quantum fields, the properties of 
which are constrained by the local structure of space and tim^. Interactions between them are 
the result of (gauge) symmetries of the Lagrangian. 

In spite of its enormous success, this cannot be a complete theory of nature for two reasons. 
On one hand, it treats gravitational fields as a classical background, while matter and other 
interactions are described by quantum field theory in the SM. This approximation certainly 
becomes invalid and has to be extended to a theory of quantum gravity at energies near the 
Planck scal^ Mp = 1.22 x 10^^ GeV. We do not address this problem in the following because 
it is of little relevance for experiments in foreseeable time. On the other hand, there are four 
experimental and observational facts which cannot be understood in the framework SM+GR. 
Three of them are widely believed to be related to particle physics, 

(I) flavour violation in neutrino experiments, section [21 

(II) the cosmological origin of the baryonic matter in the universe, section [5l 

(III) the composition and origin of the observed dark matter (DM), section [6l 

In addition to the above evidence for the existence of "new physics" , there are a number of hints 
in experimental data that may point towards the existence of physics beyond the SM; these have 
not (yet?) led to a claim of discovery and may also be explained by systematics. Of these, we 
will only discuss two in detail in this reviewH 

(i) the statistical preference for additional relativistic particles dubbed "dark radiation" (DrJI 
in fits to cosmological data, section HI 

(ii) the anomalies seen in some short baseline and reactor neutrino experiments, section [3l 



^The basic laws of other areas of natural science and technology can be understood as effective theories, which 
in principle can be derived from the SM and GR. Though there exist many complex phenomena that we do not 
understand or cannot predict in detail, this lack of predictivity is almost certainly related to the complexity of the 
system rather than a lack of understanding of its basic components, the elementary quantum fields. 

■^There are 19 free parameters in the SM; these are usually chosen as six quark masses, three mixing angles 
and one CP violating phase for the quarks, three charged lepton masses, three gauge couplings, two parameters 
in the Higgs potential and a QCD vacuum angle. The neutrinos are massless in the SM. GR adds two additional 
parameters to this barcode of nature, the Planck mass and the cosmological constant. 

■^AU known elementary particles transform under irreducible representations of the Poincare group. 

''We use natural units c — h = 1. 

^Others include the long standing issues of the muon magnetic moment (see e.g. 21 for a review), the annual 
modulation in the DAMA data Q], the excess of positrons in the cosmic radiation [i, 3I, the "forward-backward 
asymmetry" observed at the Tevatron the disputed evidence for neutrinoless double /3-decay claimed by the 
Heidelberg- Moscow experiment the cosmological lithium problem Q|, unexplained features in the galactic 7- 
ray spectrum that may be related to DM and shifts in quasar absorption lines that have been interpreted as 
signatures of a varying fine structure constant [l^] . 

^The term "dark radiation" refers to relativistic particles in the early universe with no or tiny interaction with 
the SM at temperatures T < 2 MeV. 
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All of the above phenomena may be related to right handed (RH) neutrinos with different masses. 
It is the purpose of this article to summarize how they can be connected to these hypothetical par- 
ticles and review bounds from theoretical considerations, laboratory experiments and cosmology 
on RH neutrino properties. 

To complete the list, let us add the fourth piece of evidence for physics beyond SM+GR, 
which is related to gravitation of cosmology, 

(IV) the overall geometry of the universe (isotropy, homogeneity and spatial flatness), as e.g. 
seen in the cosmic microwave background (CMB). 

We do not address (IV) here, as it is not related to RH neutrinos (or particle physics in general) 
in an obvious wajEI. An intuitive explanation is given by cosmic inflation, see section [4?n Finally, 
the observed acceleration of the universe's expansion is often included in this list. However, all 
observations can currently be explained in terms of a cosmological constant A, which is simply a 
free parameter in GR. Hence, the accelerated expansion can be accommodated in the framework 
of SM+GR. The question of the microphysical "origin" of A (and its smallness) only arises when 
the SM and GR are interpreted as low energy limits of a more general theory, including a com- 
plete description of quantum gravity. To date, (I)- (III) and (IV) are the only confirmed empirical 
proofs of physics beyond the Sivif]. 

In the remainder of this section we introduce the concept of RH neutrinos and define our 
notation. In section[2]we review how they can generate masses for the known neutrinos. In section 
[3] we summarize bounds from past laboratory experiments on RH neutrino properties, discuss the 
interpretation of the observed neutrino oscillation anomalies (ii) in terms of RH neutrinos and 
comment on possible future searches. In sections [4]l6] we discuss various cosmological constraints, 
starting with a general summary in section 14.11 The perspectives to interpret the hints for "dark 
radiation" (i) in terms of RH neutrinos and reconcile them with the oscillation anomalies (ii) are 
addressed in section 14.21 Section [5] is devoted to the idea that RH neutrinos are the origin of 
the baryonic matter in the universe (leptogenesis) and possible implications for their properties. 
Section [6] discusses RH neutrinos as DM candidates. In section [7] we address the question how 
many of these phenomena can be explained simultaneously by RH neutrinos alone. We conclude in 
section [8] and give a tabular summary of possible RH neutrino mass scales and their implications 
for known and future observations in appendix Rl 



1.1 The missing piece? 

All matter we know is composed of elementary fermions with spin ^. These can be described 
by Weyl spinors, which transform under irreducible representations of the Poncaire group, and 
combinations thereof. There are two such representations, known as "left chiral" and "right 
chiral" spinors. Remarkably, all known elementary fermions except neutrinos come in pairs of 
opposite chirality, i.e. have been observed as "left handed" (LH) and "right handed" (RH) 
particle^, see figure [H For unknown reasons the interactions of the SM are such that both can 



^It has been speculated that this point may be related to the RH neutrinos' superpartners [O], see also 
*There are various aspects of the SM that may be considered "problems" from an aesthetic viewpoint or physical 
intuition, such as the hierarchy between the electroweak and Planck scale, the strong CP problem, the factorization 
of the gauge group and the flavour structure. We do not discuss these here. We also do not discuss the issue of 
vacuum stability, which seems inconclusive at this stage due to uncertainties in the top mass 

^In this article "RH" and "LH" always refer to the chirality of the flelds and not to helicity eigenstates. 
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Figure 1: The particle content of the SM. Are we missing the right handed partner of the neutrinos? 
Picture taken from p^j . 



be combined into a Dirac spinor, see appendix |Bl Neutrinos, however, so far have only been 
observed as LH particles. One conclusion that could be derived from this is that no right chiral 
"partner" for the observed LH neutrinos exists in nature. Another possible conclusion is that we 
have not seen RH neutrinos just because their interaction with other matter is too weak. Indeed 
LH neutrinos are electrically and colour neutral; in the SM they only participate in the weak 
interaction, which does not couple to RH fields. This suggests that their RH partners are singlet 
under all gauge interaction^^. Such particles are referred to as "sterile neutrinos" . 

Let us add n RH fermions z^/j^j to the SM that are singlet under all gauge interactions and 
couple to LH neutrinos in same way as RH charged leptons couple to LH charged leptons, i.e. 
via Yukawa interactions. We will refer to these fields as RH neutrinos and to the index j that 
labels them as flavour index. Then the most general renormalizable Lagrangian in Minkowski 
space that only contains SM fields and vr reads 



^''Ii^r)- 



(1) 



Here we have suppressed flavour and isospin indices. Csm is the Lagrangian of the SM. F is 
a matrix of Yukawa couplings and Mm a Majorana mass term for the right handed neutrinos 
i^R- II = {^l^^l)^ are the left handed lepton doublets in the SM and $ is the Higgs doublet. 
$ = (e*l*)^, where e is the SU{2) antisymmetric tensor, and u'f^ = Ci^ji^, where the charge 
conjugation matrix is C = ^7270 in the Weyl representation. We chose a basis where the charged 
lepton Yukawa couplings and Mm are diagonal. Throughout this article we assume that there 
are only three "active" neutrinos that are charged under the weak interaction. If there are 
more families, then these must be heavier than mz/2, otherwise they would contribute to the 
width of the Z-boson fis'], see [l^ for some recent discussion and references. 

In the Lagrangian ([TJ the fields i^r only interact via the Yukawa couplings F. In the early 
universe, when the temperature was high enough that Higgs particles were present in the primor- 



dial plasma (T > Tew ~ 140 GeV for a Higgs mass mn ~ 125 GeV [2l|, [24]), this interaction 



^"Here we refer to the SU{3) (E) SU{2) Cgi C/(l) gauge group of the SM. They can of course be charged under some 
extended gauge group, which either acts only in a "hidden sector" or is broken at energies above the electroweak 
scale. 
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allowed z//j-particles to participate in various different scattering processes. At energies much 
below the mass of the W-boson one can in good approximation replace the Higgs field by its 
vacuum expectation value v = 174 GeV. Then ([1]) can be written as 

[{MM)ij^^iyR,j + {MmTijT^j^rA (2) 

where we defined the Dirac mass matrix nriD = Fv. Thus, at T ^ Tew the only effect of the 
Yukawa interaction is the generation of the Dirac mass term m £> , and the only way how the fields 
I'R interact with the SM is via their mixing with i/l due to m£,. 



1.2 The range of RH neutrino masses 

While in the SM the Higgs mass mn is the only dimensionful parameter (apart from the Planck 
mass), the Lagrangian ([1]) introduces n new dimensional parameters in Mm- The scale(s) asso- 
ciated with these provide a convenient way to classify different RH neutrino scenarios. Various 



embeddings of ([TJ into a bigger framework make different predictions for Mm (see e.g. 23l425l| 



for a general overview and [2&428l| for recent developments in model building), but empirically 
there are only few constraints. The following scenarios are particularly motivated, a summary in 
table form is given in appendix [Al Of course. Mm can have eigenvalues in several different mass 
ranges, so that several of these scenarios may be combined in nature. 

Mm ^ 10^ GeV - This range is motivated by embeddings of ([T]) into GUT scenarios [29|, 
such as SO(IO) unification S^, HH- SO(IO) models necessarily require the existence of fi^^ For 
Yukawa couplings F of order one, RH neutrinos with masses favoured by GUT models reproduce 
the scale of observed neutrino oscillations (I) via ([9]). In addition, typical parameter values 
allow to generate the observed baryon density in the universe (H) in CP-violating decays of RH 
neutrinos, see section [5Tl The regime Mm < 10^^ GeV is favoured bj. 

Mm ~ TeV - Theoretically this mass range is interesting because it follows from a no new 
scale principle of minimality: If Mm is near the electroweak scale, the origin of both scales 
may be related. The origin of matter (II) may be explained by CP-violating z^T^-oscillations (see 
section [5^2]) or, if two vr masses are degenerate, decays (see section [5TT|) . Neutrino masses (I) 
are explained by the seesaw mechanism. From an experimental viewpoint this mass range is 
favourable because it is accessible by high energy experiments, such as LHC. 

Mm ~ GeV - If Mm has at least two eigenvalues > 2 GeV and another eigenvalue in the keV 
range, then the observations (I)- (HI) can be described by ([1]) alone, and no other physics between 
the electroweak and Planck scales is required, see section [71 Experimentally ur with GeV masses 
may be found using high intensity experiments, see section [331 

Mm ~ keV - RH neutrinos with keV masses are promising candidates for the DM (HI) , see 
section [6l 

Mm ~ eV - RH neutrinos with eV masses can provide an explanation for the anomalies (i) 
and/or (ii), which are observed in some neutrino experiments (see section [^^ and cosmological 
data (see section [1]) . 



^^Any model that contains a U{1)b-l gauge symmetry requires this for anomaly freedom. This can be used as 
an argument for the existence of i/r: The conservation of B — L in the SM is not related to a gauge symmetry. If 
there is such symmetry, then ur must exist. 
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Mm = - For n = 3 the leptonic sector exactly resembles the quark sector without strong 
interactions. In this case neutrinos are Dirac particles Then neutrino masses are generated by 
the Higgs mechanism in precisely the same way as other fermion masses, and their smallness can 
only be assigned to very tiny Yukawa couplings. Though in principle possible, this may appear 
"unnatural" unless there is a deeper reason for i^j. Furthermore, there is no known principle 
that forbids Mm for the gauge singlet fields i'4^i. This is in contrast to quarks and charged 
leptons, for which an explicit mass term is forbidden by gauge symmetry. 



2 Neutrino oscillations 

The probably strongest motivation for the existence of i^r are neutrino oscillations, the only 
processes amongst (I)-(IV) that have been observed in the laboratory. The neutrinos i'l are 
massless in the SM^^I. In the past two decades an increasing number of neutrino experiments has 
observed neutrino flavour changes, which indicate that neutrinos are massive and oscillate, see 



e.g. |26l . l37l | for reviews with many references. The experimental results can be divided into two 
categories: the standard 3-scenario (SS) of three massive neutrinos, which we discuss here, and 
deviations from it, which we discuss in section [3l 

The interactions of neutrinos in the SM are described by the Lagrangian term 

- j=w-,-'e,w- - j=^-,^.,w- - (3) 

where g is the SU (2) coupling constant and 9w the weak mixing angle. This defines the basis of 
weak interaction eigenstates (electron, muon and tau neutrino) cj If neutrinos have a mass, then 
the mass term need not be diagonal in this basis in flavour space. In the SS, the weak interaction 
eigenstates i'L,e, ^'L,^ and u^^r are superpositions of three mass eigenstates VL,i with masses m^. 
For a given momentum, these have different energies if their masses are different, and their wave 
functions oscillate with different frequencies. Thus, the flavour decomposition of a neutrino state 
changes in time. This can explain the observed neutrino oscillations. 

2.1 The standard scenario of massive neutrinos 

There are two different ways to effectively realize the SS of three massive neutrinos with no extra 
particles. 

Majorana neutrinos - A Majorana mass term of the form 

^17Zmi,iyl + h.c. (4) 



^^For 71 7^ 3 it is in general not possible to combine all vl and vr into Dirac spinors, see appendix iBl 

^■'Some speculations on such reasons can e.g. be found in [s^ - ls^ . 

small value of Mm is, however, "natural" in the technical sense because the symmetry of £ increases in the 
limit Mm — > (there is a global U{1)b-l)- 

^®We do not consider neutrino masses as a part of the SM because we do not know what the nature (Dirac or 
Majorana) or absolute scale of the mass term is. 

^^More precisely, if one considers the general form of the interaction |^i/L,a7''C/a^ei,/3W^^ in the basis where 
charged Yukawa couplings are diagonal, then the basis of weak interaction eigenstates for vl is the one where 
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can be constructed without adding any new degrees of freedom to the SM. This term, however, 
breaks gauge invariance unless it is generated by spontaneous symmetry breaking from a gauge 



invariant term hke 38|] 



hL^f^i + h.c, (5) 

where / is some flavour matrix of dimension 1/mass. The dimension-5 operator ^ is not renor- 
malizable; in an effective field theory approach it can be understood as the low energy limit 
of renormalizable operators that is obtained after "integrating out" heavier degrees of freedom. 
These can, for example (but not necessarily!), be right handed neutrinos, cf. (llOp . Therefore a 
Majorana mass term (jj]) clearly hints towards the existence of new physics although it can be con- 
structed from SM fields only. At low energies one effectively observes the SS with only three mas- 
sive neutrinos if the energy scale related to the new physics (in case of RH neutrinos the mass Mm) 
is sufficiently high that all new particles are too heavy to be seen in neutrino experiments. The 
Majorana mass term (jlj) can be diagonalized by a transformation mu = f7jydiag(mi, m2, m3)C/J. 
In the mass base, the neutrino mixing matrix Ui, appears in the coupling to VF^ in ([3]). 



Dirac neutrinos - If neutrinos are Dirac particles, the existence of uji is directly required to 
construct the mass term 

VZmDyR + h.c. (6) 

Though this means adding new degrees of freedom to the SM at low energies, it is still a realization 
of the SS, i.e. only three massive neutrinos are observed. The can be combined with the vl 
into Dirac spinors and there are only three different masses. One could say that there are no 
new particles in the strict sense, but just additional spin states for neutrinos. We refer to both 
of these scenarios as the SS. 

A bi-unitary transformation m£) = Uyd\a.g{mi,m2i'm'i)Vp can diagonalize the mass term 
([6]), with real and positive mj. The matrix Vy is not physical and can be absorbed into a 
redefinition of the flavour vector vr. One can define a Dirac spinor -\- UyVL with a diagonal 
mass term diag(mi, m2, ms). The neutrino mixing matrix Uy then appears in the coupling of that 
Dirac spinor to in ([3]). A phenomenological prediction of Dirac neutrinos is that there is no 
neutrinoless double /3-decay. In order to be consistent with observations, the Yukawa couplings 
F have to be very small {F ~ 10~^^) compared to those of the charged leptons and quarks. 



The neutrino mixing matrix - In the basis where charged Yukawa couplings are diagonal, 
the mixing matrix Uu is identical to the Pontecorvo-Maki-Nakagawa-Sakata matrix [s^, |40] and 
can be parametrized as 

= y(23)[/^y{13)f;_^y(12)^.^g^giai/2^g^a,/2^-^) (7) 

with U±s = diag(e^*^/^, l,e^*''/^). The matrices F'^*-'^ are given by 

/ 1 \ / Ci3 Si3 \ / C12 S12 \ 

V^^^'^ = C23 S23 , , = 1 , = -Si2 C12 , (8) 

V -S23 C23 / V -S13 CIS / V 1 / 

where Cij and sij denote cos(9jj) and sin(9jj), respectively. 6jj are the neutrino mixing angles, 
ai, a2 and 5 are CP- violating phases. Many parameters of the mixing matrix Uy have been 



8 



measured in recent years. In particular, two mass square differences have been determined as 
AmgQ[ = ml — ~ 7.5 x 10~^eV^ and Am^^j^ = |m| — mf | ~ 2.4 x 10~'^eV^, the mixing angles 
are 9i2 — 34°, 623 — 39° and 613 ~ 9°; the precise best fit values differ for normal and inverted 
hierarchy (but the difference is smaller than the within the la ranges) and are given in [ll. I4ll-l43| . 



see also [4J] and 45|, |4g] for recent reviews. What remains unknown are 



the hierarchy of neutrino masses - One can distinguish between two non-equivalent setups. 
The normal hierarchy corresponds to mi < m2 < m^, with Amg^j = m^ —mf and Am^^j^ = 
I2 ^ Am^^j. The inverted hierarchy corresponds to m2 > mf > m|, with 



m5 — mt ~ m| — m: 



Am^^j = m^ — m| and Am^|.j^ = mf — m| ~ m^ — m| ^> Am^^j. 

the CP-violating phases - The Dirac phase 6 is the analogue to the CKM phase and remains 
present even for Mm = 0. Global fits to neutrino data tend to prefer 5 ~ vr [41], but are 
not conclusive. The Major ana phases ai and 02 become unphysical in the limit Mm 
because they can be absorbed into redefinitions of the fields. 

the absolute mass scale - The mass of the lightest neutrino is unknown, but the sum of 
masses is bound from above as m^ < 0.23 eV j47|] by cosmologJ^. It is also bound from 
below by the measured mass squares, rrii > 0.06 eV for normal and m^ > 0.1 eV for 
inverted hierarchy. 



2.2 Neutrino masses from RH neutrinos 

In this section we discuss how sterile neutrinos described by the Lagrangian ([1]) can give Dirac or 
Majorana mass terms to the active neutrinos. It is obvious that the Lagrangian ([2]) with Mm = 
represents Dirac neutrinos with a mass term ([6]). For Mm 7^ they generate a Majorana mass 
term ([4]). If the mass Mm is sufficiently large, RH neutrinos are heavy and neutrino oscillation 
experiments can be described by an effective Lagrangian that is obtained by integrating them 
out of 

C = CsM - i^myvl with m^, = -niDM^/m^. (9) 

If the Majorana mass is above the electroweak scale {Mm ^ v) one can already integrate out the 
UR in ([U to obtain © with / = FM~^F^, 

C = CsM + ]Tl^FMiIf^^^II. (10) 

Hence, the SS with three massive neutrinos is effectively realized if Mm is either zero or 
so much bigger than the observed neutrino masses that ([9]) can be used to describe neutrino 
oscillations. If one or more eigenvalues of Mm are not that large, then the light sterile neutrinos 
appear as new particles in neutrino experiments. To explore the full range of masses, we write 
the mass term as 

^^The bound and preferred value for ^ . mi slightly change, depending on the dataset and analysis. 
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Here we have used the identity v'^Mm^ji = i^MjjUji and the symmetry of Mm- If i^r exist, they 
always generate a mass term for ul (except for F = 0) ; this is seen as a strong motivation 
to postulate their existencj^. The relative size of Mm and niD (i.e. their eigenvalues) allows 
to distinguish different scenarios. Note that this does not say anything about the absolute scale 
of which may lie anywhere between and the electroweak scale (assuming perturbative 
couplings F). 

Mm <C niD - This is a pseudo-Dirac scenario; for n = 3 the neutrinos in this case effectively 
behave like Dirac particleCj. This case is excluded unless Mm ^ 10"^ eV, as otherwise solar 
neutrino oscillations into i^r should have been observed [64]. 

Mm ~ mo - In this case Mm and should both be of the order of the observed neutrino 
mass differences. The full mass term can be diagonalized by a (3 + n) x (3 + n) matrix U as 
DJl = i^diag(mi, 7712, m3, Ml, . . . ,Mn)h('^. This can lead to rather large mixing between and 

Mm ^ mD - This is the seesaw limit, in which the neutrino mass term is described by Q. 
This scenario is discussed in detail in the following section 12.31 The seesaw limit roughly applies 
if Mm > 1 eV. 

2.3 The seesaw mechanism 

In the limit niD ^ Mm (in terms of eigenvalues), the full (3 + n) x (3 + n) mass matrix for vi 
and vr has two distinct sets of eigenvalues, n of them (M/) are of the order of the eigenvalues of 
Mm, the remaining three (rrij) are suppressed by two powers of the active - sterile mixing matrix 

e = mDM-\ (12) 

The seesaw hierarchy 9 <^ 1 separates two distinct sets of mass eigenstates. One can block- 
diagonalize (llip by expanding in 6; on one hand one obtains the 3x3 mass matrix 

= -OMmO^ , (13) 

with eigenvalues rrij; it corresponds to the operator The known bounds on neutrino masses 
impose constraints between the values of F and Mm, see figure [2j On the other hand there is 
the n X n matrix 

Mn = MM + l{e^eMM + Ml,e^e*) (u) 

with eigenvalues M/. These matrices rui, and M^ are not diagonal and lead to flavour oscillations. 
Diagonalizing them yields the mass term 

^ (vZmf'^^vl + imM^^^vli) + h.c. (15) 

^*In this case it is not clear from the viewpoint of ^ why the fields vn should be called "neutrinos", as they 
have nothing in common with the known neutrinos except being neutral. Of course, such pariah neutrinos [which 
only interact gravitationally in ([l])] can be charged in a way that justifies this classification under an extended 
gauge group that it broken at high energies. 

^^Of course there are alternative ways to generate a neutrino mass term, see e.g. [4^ - [6^ or [23 - [25l | for a general 
overview with many references. 

^"This is not only way how Dirac fermions can arise. For instance, if the splitting between two eigenvalues of 
Mm is sufficiently small, then one can form Dirac spinor by combining different flavours, see appendix [Bl 
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with 

m^'^^g = diag(mi, m2, mg) , M^''^ = (Mi, . . . , M„) (16) 
In summary, the full matrix U that diagonahzes QJt reads 



1 - ^ee^ e 
-0t 1 _ i0t0 



VL 



Three mass states have hght masses mj ~ 0[6'^M][j] and are mainly mixings of the SU(2) charged 
fields vl, 

ul({i-lee^)uL-e,^^R^ (is) 

The remaining n mass states have masses Mj ~ 0[Mj\/] and are mainly mixings of the singlet 
fields vr, 

VR = uj, (^(1 - ^e^e*)i^R + e^ui^ . (19) 

We refer to vl as active neutrinos because they take part in the weak interactions The 
vr participate in the weak interaction ([3]), but only with an amplitude suppressed by 9; hence 
they are sterile neutrinos. The matrix Un diagonalises the sterile neutrino mass matrix Mjv0 
it can be seen as analogue to Uu. The diagonal elements of are much bigger than the off- 
diagonals and very close to the entries of Mm- Therefore one can in most cases neglect all terms 
of second order in 6 and approximate Mjy = Mm, Un = 1- However, if two eigenvalues of Mm 
are degenerate, then the Un can contain large sterile-sterile mixing angles. The active-sterile 
mixing angles are determined by the entries of the matrix 9 and always small. More precisely, 
the experimentally relevant mixing between active and sterile species is given by the matrix 
witlH 

Qai = {9Ut,U (20) 
Practically, experiments to date constrain the quantities 

Ul^Y.®aiell = Y.(^.l9li (21) 

and combinations thereof. 

This is known as the seesaw mechanism 65-6^ because increasing the eigenvalues of Mm 



pushes the masses of the sterile neutrinos up and those of active neutrinos down, just as if they 
sat on a seesa^J^. If the couplings F are of order one and Mm in the range suggested by GUT 
models, then rriy roughly reproduces the observed neutrino mass splittings. Hence, the mechanism 
provides a natural explanation for the smallness of neutrino masses. However, the scale Mm of 
the seesaw is phenomenologically almost unconstrained and may be as low as 1 eV [t^ ]. 

To make the Majorana nature of the fields explicit (and get rid of the charge conjugation 
matrix in the mass term) one can describe them in terms of Majorana spinors; we define the 
flavour vector^ 

N = vr + v''r, y = VL + vl. (22) 



^^We choose the phase convention in Un such that Mjv = f/jvdiag(Mi, . . . , Mn)!!^ with Mi real and positive. 
^^In dm) the matrix UlfO'^ = 9^ appears (rather than Q) because the Nj couple to vl.c but overlap with vf, „ 
^■^The term "seesaw" is also used for several modifications of this idea, such as the "type-II seesaw" jHol . [HB - ISa ] 
"type-Ill" seesaw [H^, "split seesaw" [g^, "radiative seesaw" 0, [zH or "inverse seesaw" [sil. IstI]. 
■^■'obviously vr = PrN and vl = -PlV, where Pr,l are chiral projectors. 
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Figure 2: A schematic illustration of the relation between F and Mm in the seesaw limit mn <C Mm- 
Individual elements of the matrices F and Mm can deviate considerably from this if there are 
cancellations in if J3)) . Plot taken from [j^j . 



Obviously, the elements "Vj are active neutrinos and the Nj sterile neutrinos. The mass and kinetic 
terms then can be combined as 

liV(i^ - M^''^)N + iv(i^ - m^s)v (23) 

Up to the normalization, this looks like the Lagrangian for a Dirac field, but one has to keep in 
mind the Majorana conditions Nj = Nf, Vj = v|. 



3 Other laboratory experiments 

The production and study of Nj particles in the laboratory is in principle possible if Mj is below 
the electroweak scale. At energies <C M/, the Nr only leave indirect traces in the laboratory. 
They manifest as higher dimensional operators [73|, such as the mass term ([9]). These can lead 
to deviations from SM predictions in different observables, such as lepton number violation or 
/3-decays. These signatures provide valuable information, but are usually not specific to RH 
neutrinos. Here we list a number of experimental setups that can constrain the properties of Nj. 
So far almost all but those in section [3T] have reported negative results, i.e. only allow to exclude 
certain parameter regions. 



3.1 Neutrino oscillation anomalies 

Accelerator experiments - Some short baseline and reactor neutrino experiments have re- 



ported deviations (ii) from the SS. A more detailed review of these results can found in [7J], 



which we follow closely here. The most prominent findings come from the LNSD experiment 
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Figure 3: Left plot: The region in the mass-mixing plane preferred by the LSND anomaly (coloured 
bands) is compared to the region allowed by MiniBooNe (coloured lines) and exclusion plots by other 
experiments as indicated in the plot; plot taken from Right plot: region allowed by the combined 
reactor and gallium anomalies, taken from ^2^] . Data from both is e.g. combined in figure 4 in \8]\] . 



, [t^, which studied transitions Ufj, — ^ and saw an excess of Ve-evenis. The similar KAR- 
MEN experiment did not see such excess [77|, but could not exclude the entire parameter space 
of the LSND anomaly [t^, mostly due to its shorter baselength. The most recent results from 
the MiniBooNE experiment [t^ from — )• i/g and — )• Uf, studies seem compatible with the 
LSND anomaly, which is in contrast to previous results. However, it has been pointed out that 
much of the sig nal comes from energies below 475 MeV, where the background evaluation is 



problematic 7J]. The ICARUS experiment [80] is also sensitive to sterile neutrinos; active-sterile 
oscillations get averaged over, but lead to an energy independent enhancement of event rates. 
Together these experiments roughly restrict the parameter space for one sterile neutrino to a 
mass splitting Am^ ~ 0.5eV^ and mixing sin^ 2^ ~ 5 x 10^'^ |7J], see figure [3l 

Reactor and gallium anomalies - There are two more anomalies that can be interpreted 



as a sign for sterile neutrinos with eV masses. One is the reactor anomaly 84]. The neutrino 
fluxes from nuclear reactors appeared to be in agreement with theory until recently. It was not 
new experimental data, but more refined theoretical calculations that lead to tension with the SS 
8^ . 84 1 . However, it should be pointed out that even in these results, about 10% of the /3-decay 
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branches are not known and can only be estimated. This involves using the ILL-experiment's 
measurement for the electron spectrum [ssl . Is^ as a reference point, hence the anomaly could 
be due to a systematic in that measurement. The other anomaly arises in the calibration of the 



GALLEX and SAGE experiments [87|,l88|. 



These anomalies have resulted in a great interest in searches for eV-mass sterile neutrinos 



Global fits to all data have e.g. been performed in 8l|, l89l . |90[|. As discussed in detail in [81[, 
the situation at this stage is not clear. The accelerator anomalies come from appearance mea- 
surements (muon to electron). If they are caused by sterile neutrino oscillations, there should 
also be a disappearance, which is not observed. This tension is present in 3 + 1 (three active 
and one eV sterile neutrino) model and remains in 3 + 2 scenarios (three active plus eV sterile). 
It reduces a bit in a 1 + 3 + 1 model (where one sterile neutrino is lighter than an active one). 
The reactor and gallium anomalies, which are due to disappearance, do not show such tension 
to other experiments. They come, however, from Vf. disappearance, which is controlled by other 
parameters than disappearance. He nce, the situation remaii is p uzzling. A detailed list of 



various proposals for future experiments 9ll494| can be found in j26l]. 



3.2 Lepton flavour violation 

The most studied consequence of the lepton number violation due to Mm is the neutrinoless 
double /3-decay discussed in the follow ing s ubsection l3.31 Mm and F also mediate flavour violation 



10' 



in the charged lepton sector |53l. I95l4l04l] . leading e.g. to muon decays — )• 67 and unitarity 



violation of the PMNS matrix [10 

In the framework of ([T]), with no other physics added to the SM, it is hard to observe these 
processes. For a generic choice of parameters the seesaw relation (jl3p implies that either the 
suppressing scale M/ is too heavy or the Yukawa couplings F are too small. A class of models 
that offers relatively good chances to observe lepton flavour violation are those where the structure 
of F and Mm is such that an approximately conserved generalized lepton number can be defined 



1071 ]. In this case the scales where lepton flavour violation and total lepton number violation 
occur can be rather different. In such models even RH neutrinos responsible for the baryon 
asymmetry of the universe may give a measurable contribution to /i — )• 67 

Models that embed ([T| into an extended framework, such as supersymmetric and grand unified 
theories, often contain interactions that violate le pton fiavour or even the total lepton number. 



This includes left-right symmetric models, see e.g. llll | and references therein. It also applies to 
various "bottom up" models that e.g. include an extended scalar sector which couples to neutri- 
nos. Hence, while non-observation of /i — )• 67 and other processes can constrain A'^/-properties, 
the observation would not necessarily be a clear sign for their existence. 

3.3 Neutrinoless double /3-decay 

If neutrin os are Majorana particles, neutrinoless double-/? decays (Oz^/3/3) are possible [sH] (see 
for recent reviews). Whether Oz^/3/3-decay occurs at observable rates depends on the 



Majorana mass matrix Mm- The lepton number violation is sufficient if at least one eigenvalue 



^^Right handed neutrinos ean also mediate lepton number violation in more exotic models. For instance, fi —>■ 
may be measurable 1 09] in models where the ur are composite objects that are hold together by strong couplings 
in a hidden sector (36l . Illd ] . 
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Figure 4: Diagram for neutrinoless double beta decay. Here we have made the lepton mimber Bow 
explicit by assigning arrows to fermion lines. The "clashing arrows" in the center of the diagram are 
allowed because neutrinos and antineutrinos are indistinguishable if they are Majorana particles. If 
some Ni are light enough, they may also be exchanged instead of Vi. The amplitude for this process 
vanishes in the limit Mm — > 0. 



of Mm is larger than the exchanged momentum ~ 100 MeV 114 - lld |. see also lod . 117]. The 
Oz^/3/3-decay can be pictured as exchange of an electron neutrino that acts as its own antiparticle 
between the nuclei, see figure [H The amplitude is given by the convolution of nuclear matrix 
elements with 



+ 0e/- 



m, 



Here p'^ 



Po 



is the exchanged momentum, with |p| ~ 100 MeV. In the regimes Mj <^ 100 



MeV and Mj ^100 MeV one can make analytical approximations and approximate the inverse 
half life time for Oi^/3/3-decays 



.-1 



I' lei , 



Af7<100MeV 



7\//>100MoV 



(24) 



nej). The 



Here G^y is the "phase space factor" (for instance Gq^ = 7. 93 x lO^^yr"^ for ''^Ge ^ 

dependence in the intermediate regime is more complicated lid ]. The Aii and Aij are nuclear 
matrix elements corresponding to light active or sterile neutrino exchange. The Aij are matrix 
elements for the contribution from heavy sterile neutrinos, sometimes called "direct contribution" 
or "contact term" (the ~ indicates the different normalization conventions for both regimes). Such 
heavy neutrinos are not really "exchanged" as propagating particles; because their mass is larger 
than the exchanged momentum, they can be integrated out and their contribution can be un- 
derstood in terms of an effective operator as (jl]). The nuclear matrix elements are a source of 
considerable uncertainty. 

Usuall y lig ht neutrinos (active or sterile) with masses <^ 100 MeV strongly dominate in the 
exchange lisl ]. However, without any extra sources of L-violation beyond ([T]), the existence of a 
mass state with Mj > 100 MeV remains a necessary condition for observable rates of Oz^/3/3-decay 
even if the direct exchange of these heavy particles does not give a significant contribution to 
T~^. When only light neutrinos are exchanged one can approximately factorize the dependencies 
on nuclear and neutrino physics and parametrize the latter in terms of the effective Majorana 
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Figure 5: Curr ent b ounds on niee (here called mpp j as a function of the smallest neutrino mass as 
summarized in "NS" refers to normal hierarchy, "IS" to inverted hierarchy. The cosmological 

limit has tightened with the publication of the Planck results ^3/- 
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me 



A/j<100MeV 



(25) 



For very light Nj one has to go beyond the seesaw approximation (j20p to calculate the active- 
sterile mixing matrix 0. 

Leaving aside some tu ned speci al cases lisl ] , the contribution from active neutrinos dominates 
in (type I) seesaw models jll(j . |119 |. If all sterile neutrino masses are below 100 MeV, they usually 
give a negative contribution, i.e. reduce the rate of Oi^/3/3-decays compared to the case where only 
active neutrinos are exchanged (which in any case is undetectably small if there is no M/ above 100 
MeV). The existence of light sterile neutrinos is strongly constrained by cosmology. In particluar, 
if sterile neutrinos compose the o bserved D M (see section [6|), then the bounds on their mixing 
imply that they do not contribute 1191 . 120]. If only the active neutrino mass states Vj contribute 
to (j25p and have an inverted mass hierarchy, then there exists a lower bound niee > 20 meV, see 
figure m For normal hierarchy there is no such bound. If also light sterile states Nj contribute 

m, iifi, [iiMS 

there is no lower bound for any hierarchy. 
So far there is no clear observation of neutrinoless double-/3 decays. The only claim of a 
detection [7l. ll24 | suggests niee = 0.32ib0.03 eV 125l |. but is disputed, as it appears to be in conflict 



''The effective Majorana mass is different from the kinetic mass (^. \Uei\ m 



2^2^1/2 
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wit h oth er observations 12d - 128| . For instance, the EXO-200 experiment finds niee < 0.14 — 0.38 
eV 

Other experiments so far only put upper bounds on the rate for this process, see figure 
[5] (an overview is also given in 129 . 130( |). 



3.4 Collider searches 

Depending on the masses Mj, there are different ways to look for A'^^-signatures in collider ex- 



periments, see 13l| for a review 



Intensity frontier - If kinematically possible, the sterile neutrinos Nj participate in all pro- 





Figure 6: Example for a contribution to Ni-decay at low energies. The Ni itself can be produced 
in meson decays. Instead of the vp-up pair there can also be charged lepton-antilepton pairs or 
quark-antiquark pairs (which hardronize) in the final state if this in kinematically allowed. 

cesses that involve active neutrinos, but with a probability that is suppressed by the sm all m i xing s 



Thi s makes it p ossible to produce them in meson decays for Mj < a few GeV 



13i - [l37| . 



^ a 

see also Q3a4l4flK1. One can distinguish two ways to look for Ni. First, they could be seen 
as missing energy in the meson decays that produce them; in two body decays this allows to 
determine their mass and (via the branching ratio) their mixing. Second, the subsequent decay 
of the Ni shown in figure [6] may also be observed (e.g. as "nothing — )• leptons" process) if one 
places a detector along the beamline. Fo r som e parameter choices, it can even be possible to 

142( 1 . Sever al experim ents of th is ty pe have be en se t 



observe both events in the same detector 



up in the past [lij - llSoj . in particular CERN PS191 Hii, H, NuTeV [uj], CHARM [lisj . 



149( 1 and WA66 150]. The bounds derived from these for n = 2 are shown in figure 



NOMAD 

[TTl They are valid under the assumption that Ni have no interactions other than those in ([1]) 
below the electroweak scale. If the mass Mj is larger than a few GeV, the Nj are too heavy to be 
produced efficiently (D or B meson decays are not possible) and it is unlikely that direct searches 
in the near future can find them. If they are, on the other hand, lighter than ~MeV, then their 
Yukawa coupling F must be very small due to (jlSp and the branching ration is very small^. 



High energy frontier - Sterile neutrinos with masses up to a few TeV in principle are within 
reach of the LHC or other future high energy experiments 13ll . Il5ll4l66l | . A promising signal in 
different scenarios are A'^z-decays that involve same sign dileptons (two leptons of the same charge) 
in the final state. However, if they only interact via the Yukawa couplings F, the branching ratio 



^'^Sterile neutrinos may also be produced in the decay of tauons [l39l - ll41l |. 

^*It has been suggested that, if Mm is generated by spontaneo us sy mmetry breaking below the electroweak scale, 



it may be possible to detect even lighter (keV) sterile neutrinos [l5l| . 
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is usually too small [162[ | ; a detection is only possible if F has a par t icula r structure that leads to 



hence allows for larger 



cancellations in the contributions of different elements to 
individual entries of F. The smallness you the Yukawa coupling in ([1]) is enforced by the seesaw 
relation ([13]) and by /x — )• 67 data flo3]. 

The perspectives are much better if the vr have interactions in addition to those in ([1]). 
For instance, it may be possible to study RH neut r inos in h igh energy collisions if the scalar 
sector is extended in comparison to the SM 15ll . Il56l. Il6ll. Il64l , in models where neutrino masses 
are generat ed v i a th e inverse seesaw mechanism fl65- 167 | or if spacetime has more than four 
dimensions 163l . ll68l ]. In left - right symmetric models 169l4l7l[ | the vr are charged under a right 
chiral SU{2) gauge symmetry, which can be broken near the TeV s ale. One promising possibility 
to test such models is to search for t he associa ted gauge bosons [172l | , but also the prop erties 
of the Nj themselves can be probed [iTiJlTsj . Cu rrent bounds are given in flZi, H, the 
right VF-boson should be heavier than about 2.5 TeV In sup e r sym metric (SUSY) theories 

the seesaw mechanism may also be studied indirectly at colliders 1781 ]: if their SUSY-partner 
(sneutrino) is the lightest sparticle, then observation s of the decays of heavier SUSY-particles 
into sterile sneutrinos can constrain their prop erties 117911 . It has been suggested that Nj can 
give a contribution to the Higgs mass, see e.g. 18CHl82l |. though this is disputed 183|. F i nally , 



LHC searches can also be promising for models that cou ple N r to Z' gauge bosons 



including those where Nj are responsible for leptogenesis 186l | 



1841 . Il85l |. 



3.5 Direct DM searches 

Sterile neutrinos are a promising DM candidate, which we discuss in detail in section [H If they 
compose all DM, then bounds from X-ray and structure formation imply that their Yuka wa 



couplings are so small that they cannot contribute to the observed neutrino oscillations [187 
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This means that, if one at the same time requires RH neutrinos to explain the two observed 
neutrino mass differences, there must be at least three of them. The small coupling also makes it 
practically impossible to directly observe these particles in collider experiments due to the tiny 
branching ratio for their production. They could at least in principle be found in direct detection 
experiments that look for interactions of DM sterile neutrinos from the inter stella,r m edium with 
atomic nuclei in the laboratory, and some experiments have been sugg ested [l89l.fl90l|. However, 
such detection would be extremely challenging (most likely impossible) due to the small mixing 
angle and the background from solar and stellar active neutrinos 191|. 



3.6 Other constraints 

Several other ways to constrain Nj properties have been suggested. If the y have eV masses as 
suggested by the oscillation anomalies (ii), they should affect /3-decays (tI, [iQ^-Esi] , which has 



not been observed. They may also leav e traces in neutrino telescopes like IceCube [196 - 1991 ] 
or detectors for dire ct DM sea rches 20d ]. The effect of keV sterile neutrinos on nuclear decays 
has been studied in a detection at the current stage seems very unlikely due to the 

background. In (203!] it was found that sterile neutrinos with TeV masses can improve the fits to 
electroweak precision data@ 

^^It has been shown that the bou nds on DM sterile neutrinos can constrain the active neutrino mass spectrum 
if additional assumptions are made [l88l ]. 
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Electroweak precision data also allows to constrain models of right handed neutrinos with electroweak scale 



masses beyond ((T}, see e.g. [20; 
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4 Thermal history of the universe 



If RH neutrinos with F ^ exist, they are necessarily produced thermally in the early 
universe. At temperatures above the electroweak scale Tew ~ 140 GeV (assuming a Higgs mass 
of ~ 125 GeV (2ll. |23|). Higgs particles are present in the plasma. This allows for Nj production 
as long as the temperature is high enough [T > Mj{T), where Mi{T) is an effective mass in 
the plasma]. The relevant processes are decays and inverse decays Nj o- ^In (or $ •(-)■ Njla, 
depending on the effective masses of the (quasi) particles in the plasma [205l4207l |) and scatterings 
(such as it o Njla)- If the masses Mtv are below the electroweak scale, iV/ are produced at 
T < Tew via active-sterile mixing [The same process can be viewed as a vl-t^r. oscillation in 
the flavour basis used in ([T])]. In addition to that, there may be other production mechanisms if 
vr have additional interactions with a hidden sector or extended Higgs sector, couple directly to 
the inflaton or are charged under a gauge symmetry that is broken above the electroweak scale. 
The presence of Ni in the plasma can have different effects in the early universe, which we will 
summarize in the following. 

All astronomical observations to date are in rather good agreement with the ACDM model of 
cosmolog}0, sometimes dubbed the "Concordance Model" or "Standard Model of Cosmology" in 
analogy to the SM. The most important cosmological parameters in the context of RH neutrinos 
are the fractions of the total energy budget of the observable universe from baryons (fi^ ~ 0.049) 
and dark matter {^dm ~ 0.265 without active neutrinos) 14 Tl]. as well as the effective number 
of neutrino species N^g in the radiation dominated epoclo. The latter can be understood as a 
measure for the expansion history of the universe. The (Hubble) rate of the universe's expansion 
is given by 

= ^Gp, (26) 



where G = Mp^ is Newton's constant and p the energy density of the universe. The contribution 
of "known neutrinos plus unknown physics" to p is usually parametrized as A^eff ^ Pu, where pi, 
is the contribution from one ultrarelativistic species and N^s is the effective number of neutrino 
species. It can be identified with the actual number of neutrino species if all neutrinos are 
effectively massless and there is no other "new physics" . It is common to parametrize 



TT 



/O7 + /Oneutrinos + [new physics effects] = /O7 + N^gpu = —T^ 



l + N,g- ( — 
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4/3' 



(27) 



where T^ and are temperature and energy density of photons and pncutrinos is the energy 
density of the SM-neutrinos. The second equality holds only for T < 0.2 MeV, i.e. after electrons 
and positrons annihilated. This is the temperature regime where Aeff is practically tested by 
observations. In the standard scenario Ngg = 3 during BBN and N^s = 3.046 at the time 
of photon decouplinglfl This assumes that each neutrino flavour has two internal degrees of 
freedom. In the SM (without neutrino masses), these are particles and antiparticles. In contrast 
to charged leptons and quarks, vl do not have a spin degree of freedom in the SM: For massless 
particles, the helicity must be equal to the chirality, hence must have left helicity. If additional 



''^A brief review of the ACDM model is given in [Jj, see e.g. [208l . l20si| for detailed introductions. 

"Radiation dominated epoch" refers to the period in the universe's history when more energy was stored in 
relativistic than in nonrelativistic degrees of freedom. 

'^'^The dev iatio n from 3 parametrises a deviation from the equilibrium distribution of neutrinos caused by 



annihilation [21 
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particles were present in the early universe, this would lead to a larger value of iVeff |21ll.[2]i]. 
Each vnj adds two neutrino degrees of freedom. If n = 3 and all of these were thermalized, 
would be 6. However, if neutrinos are Dirac fields, then the Yukawa couplings F must be 
tiny to explain the smallness of the neutrino masses. Then the degrees of freedom do not get 
populated significantly by thermal production. In other words, the helicity changing processes 
are so suppressed at T ^ that practically all neutrinos are produced by the weak interaction 
and have left helicity. This argument of course assumes that have no other couplings than 
F that could contribute to thermal production. Hence, the fact that we observe N^s ~ 3 — 4 
strongly constrains models of Dirac neutrinos in which are charged under some gauge group 
at high energies or otherwise produced in the early universe. If neutrinos are Majorana fields 
and Mm ^ 100 MeV, then the heavy particles Ni have decoupled and decayed long before big 
bang nucleosynthesis (BBN) and do not affect light element abundances or the CMB. Lighter 
sterile neutrinos can be long lived enough to contribute to A'^efT during BBN and afterwards. If 
they are in thermal equilibrium and ultrarelativistic, then each Nj increases A'^eff by one. If their 
abundance is below equilibrium or their mass not negligible, they contribute less. 

Since (j27p assumes that all neutrinos are massless and in thermal and chemical equilibrium, 
any deviation of the momentum distribution from {e^/'^ + [such as a chemical potential, non- 
negligible mas^fll or a nonequilibrium distribution] leads to non- integer contribution to A'eff ; this 
applies to both, active and sterile neutrinos. Furthermore, in the SM the ratio of the temperatures 
of the neutrino and photon backgrounds at the time of photon decoupling is (4/11)^/^. If this 
ratio is different (e.g. because some decaying particle injects energy into either of t hem) , the n 



this would lead to deviations from A'eff = 3.046 even if there are only three neutrinos [216l4218l |. 



4.1 A brief history of the universe 

Observations of the CMB show that the universe was homogeneous and isotropic to one p art in 



~ 100000 at redshift z ^ 1100, when photons decoupled from the primordial plasma [2191 ]. This 
is puzzling, as the radiation we receive from different directions originates from regions that were 
causally disconnected at that time if the universe only contained radiation and matt er (" horizon 



problem"). Furthermore, the inferred overall spatial curvature is zero or very small liia ]. which 



means that it was extremely close to zero at earlier tim e s ("f latness problem"). Both problems 



can be understood as the result of cosmic inflation [220l . l221[, a phase of accelerated expansion 
in the universe's very early history. Inflation can also explain the small density perturbations 
that served as seeds for structure formation in the universe as quantum fluctuations that were 
"streched out" to macroscopic scales by the rapid expansion, and predict the correct properties 
of their spectrum. However, while inflation is an excellent model for cosmology, we do not know 
much about the fundamental physics mechanism that drove it. 

If inflation was driven by the potential energy of an inflaton field, then the quantum fluc- 
tuations of this field lead to small perturbations in the (otherwise homogeneous) gravitational 
potential. After the inflaton dissipates its energy into relativistic particles during cosmic re- 



heating these lead to density fluctuations in the primordial plasma, which manifest 



in temperatu re flu ctuations in the CMB and form the seeds for the formation of structures in 
the universe The power spectrum of CMB fluctuations is in very good agreement with 



See [2131 . |214| | for a detailed discussion of the effects of neutrino masses in cosmology. For eV masses the 
deviatio n is small in the early universe. An early discussion in Russian can be found in [2151 ]. 
^^See [22511 for a recent discussion. 
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the above hypothesis 219l |. After reheating, Nj can affect the thermal history of the expanding 
universe in several ways. 



Baryogenesis via leptogenesis (T > Tew) - There is good evidence [227l | that the observed 
is the thermal rehc of a smah matter-antimatter asymmetry in order ~ 10"^'' of the early 
universe (BAU), which survived after all other particles and antiparticles annihilated into CMB 
photons and neutrinos and is reflected in today's baryon to photon ratij^ rjB ~ 10^^'^. If this 
asymmetry was produced by leptogenesis, then the observed CIb allows to constrain F and Mm, 
cf. section [5l 



Dark Matter production (T ~ 100 MeV if produced by mixing) - If DM consists of thermally 
produced Nj, then there are various different bounds on their properties, which we discuss in 
section [6l 



Late time leptogenesis (Tew > T > few MeV) - Sphaleron processes are the only 

source of baryon number violation in ([T]). They freeze oull^ at T ~ Tew-, below which baryon 
number is conserved. The production of a lepton asymmetry can, however, continue afterwards 
if some Nj are out of equilibrium (e.g. during their freezeout and decay) because F and Mm 
violate flavoured and total lepton number, respectiv ely. The generated asymmetries can be much 
bigger than ijb and differ in each flavour 230, l23ll |. A weak cons traint on the asymmetry may 



be derived from its effect on hadronisation at T" ~ 200 MeV |23i |. Stronger constraints can be 
derived from BEN, see below. Late time asymmetries can be very important if there are long 
lived sterile neutrinos because they can amplify or suppres s th e ir pr oduction rate. For keV-mass 
sterile neutrinos an amplification of their p rodu ction rate 233l . 12341 ] can be so efficient that they 
are abundant enough to constitute all DM 230l |. A suppression of the thermal production rate, 
on the other hand, can help to ease the tension between cosmological and laboratory hints for 
eV-mass sterile neutrinos, s ee se ction 14. 2[ In return, light sterile neutrinos can also amplify an 
active neutrino asymmetry [235| |. 



Neutrino freezeout (T ~ 1.1 MeV for active neutrinos) - The freezeout of active neutrinos 
leads to a cosmic background of relativistic neutrinos, analogue to the CMB. They affect the 
expansion by their energy density /Oneutrinos- The neutrino background may carry a lepton asym- 
metry that is orders of magnitude larger than t he b aryon asymmetry ~ 77^ ~ lO^^*'. The main 
constrains on such asymmetry come from BBN [236[ . In the SS, a ctive neutrino oscillations tend 



to make the asymmetries in individual flavours equalffi 237h239( |. How complete this equilibra- 
tion is depen ds on the mixing angle 613. The measured value for 613 suggests a high degree of 



equilibration 24d |. Lepton asymmetries in the neutrino background can be constrained due to 
their effect on the momentum distribution, which changes the relation between temperature and 



The relation between these parameters is given by rjB ^ 2.739 ■ IQ-'^h^nB- 

■^^A reaction "freezes out" when the temperature dependent rate at which it occurs falls below the Hubble rate 
H of cosmic expansion. This essentially means that the density of the primordial plasma has become so low that 
the rate at which the participating particles meet is negligible. A particle freezes out when the last process (other 
than decay) that changes its comoving number density freezes out. 

■^^This is often referred to as "flavour equilibration", though the process occurs close to the neutrino freezeout 
and neutrinos may not reach thermal equilibrium. This means that the lepton asymmetry cannot be translated 
into a chemical potential in the strict sense. 
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energy density ()27p and mimics an A'cfr 7^ 3. T he b ounds on N^g from BBN (see below) allow to 
cons train the asymmetry to roughly \r]L\ < 0.1 [2401 ] . where tjl is defined analogously to ([281) : see 
also |24ll - [243 |. 

If sterile neutrinos are relativistic, then they can form a similar background and contribute 
to A'efr. For an active-sterile mass splitting Am^ < 1.3 x 10~^eV^ active-sterile oscillations 
are effective only after the active neutrino freezeout, then they simply distort the momentum 
distributions 214[ |. For much larger splittings they can be produced efficiently via their mixing 
at T > 1 MeV. The preference for A^eff > 3 in different cosmological data sets (see below) 
can be interpreted as a hint for eV mass sterile neutrinos as DR, see section 14. 2[ However, the 
constraints on from BBN, see below, imply that such background must either have frozen out 
considerably earlier if it ever was in equilibrium or never thermalized (e.g. because the production 
was suppressed by some mechanism). 



Big bang nucleosynthesis (T < 100 keV) - There was a brief period in the early universe 
during which the temperature was low enough for nuclei heavier than hydrogen (H) to exist and 
still high enough for thermonuclear reactions to occur. During this period most of the deuterium 
(D), helium ('^He, ^He) and lithium (mainly ''Li) in the universe were formed [245i ]. see e.g. chapter 
22 in [l[ for a review. These light elements, in particular ^He, make up the vast majority of all 
nuclei other than H in the universe. Sterile neutrinos can affect this big bang nucleosynthesis 
(BBN) in different ways, depending on their mass. Nj with masses far above the electroweak scale 
have no effect on BBN, as they have decayed long before. If the masses are in the GeV to TeV 
range, these particles can be long lived enough that the entropy released during their decay affects 
BBN or the thermal history afterwards. The good agreement between BBN calculations and the 
observed H and He abundances implies that, if sterile neutrinos with Ge y< Mr ^-.TeV exist, 
they must have decayed sufficiently long before BBN. The r esulting bounds |23, EIEi] in the 
mass-mixing plane are plotted in figure [TTJ If DM is composed of keV-mass sterile neutrinos, 
these would have no visible effect on BBN because their number density is too low to affect the 
expansion history in that era. Li ght s teril e neu trinos (with eV masses) would, on the other hand, 
significantly contribute to N^f^ 213, 247 . 248| as additional radiation and increase the rate of 
expansion of the universe via (I26p and (I27p . This, for instance, determines the precise moment 
of neutron freeze-out (T ~ 0.8), which roughly occurs when the expansion rate (|26p equals the 
rate for the reaction n + e"*" o p + z^e ■ It also affects the amount of time that passes until the 
formation of elements (T ~ 10 — 100 keV) and duri ng which neutrons decay. Both determine the 



number of neutrons available for fusions 249], |250|. However, a change in the expansion rate is 



not specific to light neutrinos, hence the resulting change in the value for A'efr extracted from fits 
to light element a bund ances does not actually " measure" the number of light neutrinos or other 
particles (see e.g. 25l[ | for a summary of examples). It simply parametrizes any deviation from 



the SM prediction. Sterile neutrinos can also affect BBN if they cause deviations of the neutrino 
momentum distributions frorna Fermi-Dirac spectrum, e.g. by active-sterile oscillations or by 
inducing chemical potentials 253] • On one hand, this slightly modifies the relation ([27p between 

T and p,^, which again affects the rate of expansion.^ More importantly, the He abundance is 

directly affected by a distortion of the sp ectrum 213| . D is less affected, hence provides a 
more direct probe of the expansion rate 2531 ]. 

BBN depends on Q^, N^q and the lepton asymmetries Y^. In the SM, where A'efT = 3 is fixed 
and Ya = 0, the only free parameter during BBN is (or tjb), and there is a rather impressive 
agreement between theoretical predictions and the value obtained from the light elements, see 
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e.g. 

in the 
also 



254l |. If one treats N^q as a free parameter using (j27|) . a deviation from 3 mainly refl ects 
^H e ab undance. The values for N^g obtained from BBN alone (see e.g. 254 - 256| . cf. 
257 . 258| |) show a slight preference f or N pff > 3, but are consistent with A'efr 



3. Different 

interpretation of the data are discussed in 2501 ] : In absence of significant lepton asymmetry, BBN 



3.71 



+0.47 



If one fixes the He abundance to the value inferred from the CMB. 



alone yields N^s — ^-i ^~0 45- 
this tightens to A^cff = 3.53lo;g3 (A'efi = 3.22 ifc 0-55 when using the D abundance measured in 
2561 ] alone). The likelihood functions given in [255] show a preference for Nes > 0, but clearly 
constrain N^s < 4. Hence, sterile neutrinos with eV masses, as motivated by the oscillation 
anomalies (ii), can only be made consistent with BBN if they are not thermalized. This strongly 
constrains models of light sterile neutrinos unless some mechanism suppresses their production. 
If DM is composed of sterile neutrinos with keV masses (see section [6]) , these do not affect 
because their number density is well below the equilibrium value (their number density can simply 
be obtained by putting today's 17^ in relation to their mass). Hence, they are unaffected by these 
bounds. 



Matter - radiation equality (T ~ 0.8 eV) - As the universe expands, the matter density is 
diluted as oc a'^, where a is the scale factor, due to the increasing physical volume. The radiation 
energy density is diluted faster (oc a"^) due to the stretching of the wavelengths, hence an initially 
radiation dominated universe becomes matter dominated at the point of matter-radiation equality. 
In cosmology, one refers to all relativistic degrees of freedom (particle energy dominated by 
momentum) as "radiation" , while nonrelativistic degrees of freedom (particle energy dominated 
by mass) are "matter". A first principles derivation of the kinetic equation in the expanding 



universe that covers both regimes is given in 259|]. Due to Hubble expansion, constituents 



of the primordial plasma change their identity from "radiation" to "matter" when T (more 
precisely: their average momentum) falls below their mass, which has to be taken into account 
when determining the point of equality. If sterile neutrinos are relativistic near T ~ 1 eV, they 
change the temperature of matter radiation equality. This is crucial for the growth of density 
perturbations in the primordial plasma, see following paragraph. If they become nonrelativistic 
just around this time, their equation of state is neither that of "radiation" nor "matter" in the 
intermedia te re gime. The effect of the time dependence of the equation of state has e.g. been 



studied in [2601]. 



Photon decoupling (T ~ 0.25 eV) - The universe becomes transparent when the temperature 
is so low that photons cannot dissociate H atoms any mor^ff]. While photons previously scattered 
frequently with free electrons, the cross section with neutral atoms is so small that the average 
photon has not interacted with matter ever since. The cosmic microwave background (CMB) 
formed by these primordial photons allows us to observe the universe at a very early stage 
(redshift z ~ 1100). The CMB contains an enormous amount of information and is, altogether, 
one of the most impressive confirmations of the ACDM model. The vast amount of data also 
makes it an excellent tool to look for hints of physics beyond ACDM and the SM. The small 
temperature fiuctuations (of relative order 6T/T ~ 10~^) in the CMB yield the earliest probe of 
structure in the universe. They were generated by acoustic oscillations of the coupled baryon- 

''^The temperature at which this happens is much lower than the H binding energy because of the small baryon 
to photon ratio rjB, which implies that there are enough photons in the high energy tail of the Bose-Einstein 
distribution to dissociate the atoms. 
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photon plasma in the primordial gravitational potential wells (DM also falls into these wells, but 
does not feel the radiation pressure and does not oscillate). 

The fluctuations can be decomposed into modes with wave numbers k, their spatial extension 
is characterized by the inverse of k. When is larger than the causal horizon, the mode remains 
"frozen". After the end of inflation, density fluctuations successively "(re)enter the horizon"; that 
is, the causal horizon becomes bigger than their spatial extension. Then their amplitude starts to 
grow due to gravitational infall. The evolution in the radiation dominated era is governed by the 
competition between gravity and the radiation pressure (which depends on Nes); this competition 
leads to haryon acoustic oscillations (BAO) of the plasma. The modes oscillate under the action 
of both forces until they decouple (when the baryons stop feeling the pressure) and start to 
collapse under the action of gravity, eventually forming the the structures that we observe in the 
universe. The observed peaks in the CMB power spectrum can be identified with multipoles / in 
the multipole expansion of the CMB fluctuations that correspond to the modes k which reached 
the maximal elongation at the moment of decoupling. The first peak corresponds to the smallest 
k that reached maximal compression by the time of decoupling and so on. The evolution of 
perturbations can be studied quantitatively by a coupled set of Boltzmann and Einstein equations. 
Before photon decoupling, the density perturbations are small and 5kp/ p ^ S can be used as 
as an expansion parameter in calculations. 

T he powe r spe ctrum of temperature fluctuations can be affected by RH neutrinos in different 



ways [261. |214| . l26l[ |. If they are nonrelativistic, they act as DM. Then they only affect the growth 
of perturbations via their free streaming length (the DM sterile neutrinos discussed in section 
[6] are relativistic at freezeout, but become nonrelativistic between BBN and the decoupling of 
photons). If they are relativistic (as e.g. the eV-mass-A^/ motivat ed by os cillation anomalies). 



they act as radiation and modify the power spectrum in several ways 26ll4267l |. Their contribution 
to Aeff > 3 increases the rate of expansion ()26p . hence reduces the comoving sound horizon, which 
moves the peaks in the power spectrum to higher m ultip oles. It also enhances the hight of the 
first two peaks via the integrated Sachs- Wolfe effect |263| . Light sterile neutrinos are essentially 
collisionless and can, in contrast to baryons and photons, not be treated as a perfect fluid. Their 
anisotropic stress affects the gravitational potential via Einstein's equations. Since this effect is 
more relevant in the radiation dominated era, it affects modes that enter the horizon before and 
after matter-radiation eq ualit y in a different way, which reflects in a change in the relative height 



of high and low / peaks [2651]; it suppresses the modes / > 200 and also shifts the positions of 
the peaks [26 li] . Finally, sterile neutrin os cq r itribute t o the damping of small scales via their free 
streaming and affect the high-/ modes 261 . 262, 26^. All these effects can be parametrized in 



terms of A^effQ 

The CMB power spectrum has been studied with different instruments, including the Wilkin- 
son Microwave Anisotropy Probe (WMAP) satellite, South Pole Telescope (SPT), Atacama Cos- 
mology Telescope (ACT) and Planck satellite, see figure [71 A few years ago these datasets 
consistently preferred Agfj > 3. In order to constrain Aeff , it is crucial to have information from 
high and low multipoles /. While WMAP with its full sky coverage can measure the low I at 
good precision, its angular resolution is not high enough to go beyond the third peak. ACT and 
SPT have a better resolution, but can only observe parts of the sky and cannot go below / ^ 500. 



'"^Here &pk/ p collectively refers to the contrasts in the densities of photons, baryons, neutrinos and metric degrees 
of freedom, which of course have to be studied independently. 

*^Note that individually they are degenerate with other effects, cf. e.g. the discussion in and references 
therein, and only a global fit allows to reliably determine A'eff . 
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Figure 7: CMB power sgectrum as measured by WMAP 12191] . SPT [269] . ACT \27d] and Planck 
plot taken from pT^] . 



Thus, these data sets need to be combined to reduce the degeneracies between N^.^ and other 
parameters. Planck has, for the first time, measured low and high I with a single instrument. 
Combining WMAP 7 year data with ACT in 2011 preferred A'efr = 5.3 it 1.3 (which reduced to 
4.6 ±0.8 if data from BAO and measurements of the Hubble rate Hq are adde d) f27oll . combining 
SPT with WMAP7 gave iVeff = 3.85 ± 0.62 (3.86 ± 0.42 with BAO and Ho) |26gt]. More recent 
SPT results [272t] continue to favour DR with iVcff = 3.62 ±0.48 (3.71 ±0.43 when combined with 
BAO and Hq data). In contrast to that, new ACT data [ili, favours iVeff = 2.79 ± 0.56, in 
accord with the SS. This tension between SPT and ACT has e.g. been discussed in 0,[273. The 
WMAP colla borat ion quotes N^^ = 3.84 ±0.40 using their 9 year data |2l3- A combined analysis 



performed in [276|] yields TVeff = 3.28 ± 0.40. Very recently, Planck found N^s = 3.30 ± 0.27 [43], 
which is perfectly consistent with the SS. We discuss the implications of these measurements for 
sterile neutrinos as DR in section 14.21 



Large Scale Structure The formation of structures in the universe (see 278| for a review) 
after photon decoupling can roughly be divided into two parts. As long as 5pk/p <C 1 can be used 
as an expansion paramet er (" l inear regime"), seniianalytic methods [279] and fully automated 



codes (such as CAMb|11 [isij or CLASdl!! isJl) are applicable. These can test a large 



''^ http : / /camb . info 

■^"^ http://lesgourg.web.cern.ch/lesgourg/class.php 
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Figure 8: Matter power spectrum from diSerent observations as indicated in the plot and combined 
in ',277,] (plot taken from there). 



number of cosmological parameter sets in a short t ime. For 5pk/p > 1 ("nonlinear regime") 
the apphcabihty of semianalytic metho ds (see e .g. 2831 - 285 |) is very hmited; in general only 
expensive numerical n-body simulations 28a- 292] allow to make predictions. In the intermediate 
range, some sophisticated methods have been suggested (see e.g. [293,- ,299l] ) . For perturbations 
on different scales, the evolution becomes nonlinear at different times; the behaviour on cosmic 
scales is still linear nowadays, while the matter distribution locally is extremely inhomogeneous. 
Sterile neutrinos can affect the structure formation if they act as DM or DR. 

The implications of DM sterile neutrinos with keV masses are discussed in section [6l the 
only big difference to CDM scenarios lies in their free streaming length Xdm- This leads to 
a suppression of structures on (comoving) scales smaller than Xdm- Lighter sterile neutrinos 
with ~eV masses as suggested by the oscillation anomalies (ii) would act as DR and as HDM 
contribution to ^dm- The effect that they have on the expansion rate ()26p would not only reflect 
in the abundances of light elements and the CMB, but also in the distribution of matter in the 
universe. 

An important quantity that characterizes the distribution of large scale structures (LSS) in the 
universe is the matter power spectrum, cf. figure El which is obtained from large galaxy surveys 
and measures fluctuations on smaller scales than the CMB. Apart from small wiggles 



(related to baryonic acoustic oscillations), the matter power spectrum shows a "turnover point", 
the position of which is sensitive to N^s- This feature occurs because perturbations oscillated in 
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Fi gure 9: Left plot: Combinations of Ncs and the sum of neutrino masses favoured by the analysis 
in |9d /. Note that the parameter Ncs on the horizontal axis is not exactly identical to the parameter 
in i|271) because the authors allowed for a common neutrino mass. Right plot: Marginalized posterior 
distribut ion of NcH from Planck (black) a nd P lanck data supplemented by BAO data '303 , 305 , 30a. 
\30A . \3ldij (blue), a local Ho measurement \3ld l (red), and both BAO and Ho (green); plot from [47l /. 



the radiation dominated era after they enter the horizon and grew only slowly, while in the matter 
dominated era they grew more quickly. The initial power spectrum is almost scale invariant, but 
those modes that entered the horizon during radiation domination grew only slowly until the 
moment of matter-radiation equality and are suppressed with respect to those that enter after 
that point. A larger A'efr leads to a later equality, which affects the position of the turnover point. 
The damping of modes right of the turnover point is controlled by the ratio Qs/i^B + ^dm)- It 
is also affected by the velocity dispersion of massive neutrinos [2g| . However, since the position 



of the turnover point is only known with limited accuracy 31l[ and LSS data cannot constrain 



all cosmological parameters, it is usually analysed in combination with CMB data (see below). 
4.2 Sterile neutrinos as dark radiation 

For the past decade, cosmological data generally favoured N^s > 3.046, though the significance 
was not strong enough to rule out the SS. It is tempting to think that sterile neutrinos may be 
the common origin of the oscillation anomalies (ii) and these observations (i) . For the masses and 
mixings indicated by LSND, the thermal pro duction of the sterile neutrinos in the early universe 



would be efficient enough to thermalize them 312l4314i |. leading to a contribution to iVefj of order 
one per species. 

Unfortunately, the current situation is rather confusing. On one hand there is some tension 



amongst different oscillation experiments 81[; most of the region favoured by LNSD is excluded 
by ICARUS, KARMEN and MiniBooNE. On the other hand, there is also tension between the 
latest results from the CMB. SPT+WMAP results [272,] favour DR with N^s = 3.62 ± 0.48; 
this is consistent with the A'^cfr — 4 predicted by a thermalized sterile neutrino with masses and 
mixing suggested by the oscillation anomalies, but also does not rule out the SS. In contrast. 
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ACT+WMAP Hzi, Szil favours N^s = 2.79 ± 0.56, which would clearly rule out thermal DR. 
This is in contrast to earlier ACT results 27Cll |. which favoured DR. Planck found 3.30 it 0.27 
[47! (also using BAO data); though this shows a slight preference for A'eff > 3, it is perfectly 
consistent with the SS and clearly disfavours A'efT = 4. 

The WMAP, ACT and SPT results have been combined with various other data sets, including 
measurements of the CMB, the LSS and BAO from galaxy surveys and measurements of the 
contemporary Hubble constant Hq. For instance, in j90t| it was recently argued that cosmology 
can be consistent with the neutrino oscillation anomalies. This conclusion was mainly possible 
because the authors excluded ACT and BBN data from the combined analysis. This can be seen 
in figure [H For small masses < 0.5 eV, WMAP+SPT+ACT data prefer extra neutrinos; but 
for masses around leV the data favours A'efr ~ 3, as predicted by the SS. Including BAO and 
Hq data indicates N^fi > 3, but also disfavours the ~ 1 eV region for the masses preferred by 
oscillation anomalies. Without ACT, = 4 and a mass of 1 eV are perfe c tly consis tent. Similar 
analyses, using pre-Planck data, have e.g. been performed in |258l. [iei. \27di. [316l-[322pl. The 
conclusions differ, depending on the selected data sets and statistical method. There is a general 
preference for A^gg > 3, but no exclusion of the SS. 

However, while the authors of [90 ] presented arguments to exclude ACT from their analysis, 
the conclusions change dramatically when the Planck result 3.30±0.27 is included. Though a 
slight preference for Aefj so far seems to remain, the thermalized DR {Ncs > 4) that the oscillation 
anomalies hint at is clearly disfavoured The above result assumes massless neutrinos, but 

does not change significantly if the sum of neutrino masses is allowed to vary {Nes = 3.32^q52 
and ^jiTii < 0.28 eV when combined with BAO 47|]). There is some freedom to increase A^eff if 
local measurements are used to determine Hq instead of the value obtained from the Planck data 
itself, see figure [H Local measurements of cosmological quantities such as the Hubble rate or the 
age of the universe (by studying old objects) depend les s on the assumed model of cosmology and 
can be used complementary to high redshift data, see for a recent discussion. The Planck 
collaboration concludes that the tension between direct Hq measurements and the CMB and 
BAO data in the ACDM model can be relieved by N^s > 3.046, but there is no strong preference 
for this extension from the CMB damping tail |47|. It remains to be seen whether this tendency 
continues when the full Planck d ata is analysed, which will provide the most stringent constraint 
on Nf.fl from a single instrument 3261 ]. 

The apparent tension between oscillation anomalies and the Planck and ACT results could 
be eased if the thermal production of ligh t ste r ile neutr inos in the early universe is somehow 

suppresse d, e. g. b y a chemical potential ^Ei, Bjsigl ] . Quantitative studies of the flavour 

evolution ^M, [237p239l . BB ] are required to answer the question to whic h degree a lepton 
chemical potential can prevent the thermalization 258l . |314| . l328l . l329l . 133314341 ] . 

If one, on the other hand, interprets the Planck and ACT results as indicators that there is no 
DR, then it remains to be understood why so many previous analyses seemed to point towards 
its existence. The possibility that the preference N^s > 3.0 46 is an artefact of the prior choice 
in B ayesian analyses has been discussed and rejected in 342| . On the other hand, the authors of 



3221 ] found that Bayesian model selection prefers SS, despite the fact that parameter estimates 



f^ATeff are larger than 3.046. In [343] it has been pointed out that there is no evidence for 
Acff > 3 if one relaxes the assumptions on the late history of the universe suggested by ACDM 
when analysing the CMB. Generally, there is a lot more freedom if one allows for more physics 



A nice summary of various other ways to combine different data sets can be found in 1323 
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beyond the ACDM scenario. In [217i | it has been discussed in detail that different modifications of 
the SM can predict values N^s > 3 or A'efr < 3, irrespective of the number of neutrino degrees of 
freedom. However, it seems that at least the simplest scenario in which light sterile neutrinos with 
generic couplings are the common explanation for the oscillation anomalies and the preference 
for iVgff > 3 is disfavoured by Planck unless more new physics is added. At the same time, it 
seems unlikely that the SS can be ruled out by CMB observations in foreseeable time. 



4.3 Sterile neutrinos in astrophysics 

In addition to possible signatures in high redshift observations, sterile neutrinos may also have an 
effect on astrophysical phenomena at present time. Their most studied role in astrophysics is that 
of a DM candidate, see section 14.21 but they may have effects on other phen omena . They c an, 
for instance, affect the transport in supernova explosions if they have eV 344], keV 3451 . 346 1 or 
GeV 3471 ] masses. Sterile neu trinos with keV masses can also help to explain the high rotation 
velocities of pulsars [348l . 13491 ] . 



5 Baryogenesis via leptogenesis 



The observable universe does not contain any significant amounts of antibaryons 2271 ] . i.e. it 
is highly matter-antimatter asymmetric. Given our knowledge about the thermal history of 
the universe, today's baryon density is easily explained as the remnant of a small matter- 
antimatter asymmetry at early times, when the temperature was high enough for pair creation 
processes to occur faster than the Hubble rate. The baryon asymmetry in the early universe 
(BAU) can be estimated by the baryon-to-photon ratio at later times. 
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Here ub and are comoving numb er de nsities for baryons and anti-baryons. rjB can be deter- 
mined rather consistently from BEN 254 ] or the CMB and LSS 219l ]. 
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(29) 



both in units of 10~ . A period of cosmic inflation, as suggested by the CMB, would have 
diluted any pre-inflationary asymmetry. Thus, the BAU must have been created dynamically 
afterwards. There are three con ditio ns for the dynamical generation of a BAU ("baryogenesis"), 
known as Sakharov conditions [350*]: Baryon number (B) violation, breaking of charge (C) as 
well as charge-parity (CP) symmetry and a deviation from thermal equilibrium. In prin c iple, all 



3511 ] at 



prmc 

of them are fulfilled in the SM: at the quantum leve l, baryon number is violated |228l . 
T > Tew by rionperturbative sphaleron proc esses 122911 . C and CP are violated by the phase in the 
CKM matrix [352] and the weak interaction [3531 . l354i ] , and a deviation from thermal equilibrium 



is caused by Hubble expansion. However, the CP violation and deviation from equilibrium are 
both too small in the SM to explain the observed r]B', see [2271 ] for a detailed account of the 
Sakharov conditions and baryogenesis in the SM. 

RH neutrinos described by ([T|) can fix both of these shortcomings. As gauge singlets, they 
can be out of equilibrium at temperatures when all other particles are tightly coupled by gauge 
interactions, and their Yukawa couplings contain several unconstrained CP-violating phases. This 
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makes baryogenesis via leptogenesis possible [3551) I . Practically leptogenesis works as follows. 



As usual, we define the left handed lepton numbers Li^a as the zero components of the currents 
= T^L,al^i'L,a + (iL,a7^eL,a- Analogously, we define lepton numbers Lrj and L_R,a from 
currents Jj = i'rj'J^i^rj and Jji^ ~ ^R,al^^R,a for the RH fields. This allows to define the 
active lepton numbers = Ll^o + Lr^ a well as the LH and RH lepton numbers Lj^ — ^L,a 
and Lr = ^jLrj + '^a^R,a- '^^^ total lepton number is L = + Lr. The interactions in 
([1]) violate individual lepton flavours La (via F) as well as total lepton number L (via MmJIEI. 
When one or more of the Nj are out of equilibrium, they can generate lepton asymmetries via 
various processes, including decays and inverse decays, scatterings and flavour oscillations (the 
probability to decay/scatter /oscillate into leptons and antileptons is different at second order in 
F due to quantum interferences between CP-violating processes). Sphalerons violate B and L 
individually, but conserve B — L. At T > Tew they transform part of the lepton asymmetr y int o 
a baryon asymmetry; if all interactions are in equilibrium it is given by S ~ ||(S - L) [353, l358|. 



B is protected from washout after sphaleron freezeout around T ~ Tew- Since sphalerons only 
couple to left chiral fields, they in fact only see the leptonic charge Ll stored in these fields, which 
may differ from the total lepton asymmetry L. Thus, the BAU we observe is determined by L^, 
at T ~ Tew ~ 140 GeV. The nonequilibrium condition can be fulfilled three times for each Nf. 
during their production, their freezeout, and their decay. 

5.1 Leptogenesis from Nj freezeout and decay 




Figure 10: A crucial contribution to the CP-violation in the simplest version of "vanilla leptogenesis" 
comes from the quantum interference between the tree level diagram and radiative corrections to the 
decay Ni — >■ $Zq, . Due to the CP-violation contained in the complex phases in F , the amplitudes 
for these processes with la in the final state differ from those with la. At tree level this difference 
amounts to an overall phase, which does not affect the physical absolute square of the amplitude, but 
in the interference terms between tree level and loop diagrams it makes a physical difference. The 
rightmost diagram gets resonantly enhanced for Mi ~ M,j because then the Nj propagator is nearly 
onshell. In more general scenarios, also scattering processes (such as ^la — >■ with intermediate 
Ni ) have to be considered. 

It is common to split the lepton number violating processes into "source" and "washout". 
The source consists of all contributions that violate active lepton numbers in the absence of 
an existing asymmetry (e.g. A'^/-decays) ; these can generate asymmetries in the presence of a 
deviation from thermal equilibrium. The "washout" consists of all other processes, which tend 

■^^The term leptogenesis is often used in a wider sense, relerring to all scenarios where the source ol CP-violation 
lies in the leptonic sector. Here we only refer to those scenarios where the BAU is generated from the interactions 
in IDl. 

''^If F and Aim have some particular structure, a combination of them may be exactly or approximately conserved, 
in which case a different generalization of the term "lepton number" to the RH neutrinos may be convenient, see 
e.g. tlOS, 356]. 
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to eliminat e ex i sting asymmetries (e.g. inverse A^/-decays) . In addition, there are "spectator 



processes" [359|, 1360] that do not violate L themselves, but redistribute charges amongst different 



fields and thereby affect the time evolution of L. If the Nj have masses as suggested by GUT 
models, they are produced thermally, freeze out and decay at temperatures T 3> Tew- Since 
L- violating processes are strongly suppressed at T < M/, B — L ^ Q \s preserved from washout. 
This is the most studie d version of leptogenesis, and there exist several detailed reviews, see e.g. 



36l| or, more recently, [362l . l363l | . This setup is very appealing because it can easily be integrated 
into a GUT-framework, provides a "natural" explanation for the smallness of the neutrino masses 
rrii via the seesaw mechanism and allows to probe at least some parameters of very high energy 
physics in low energy neutrino experiments via (j9]). The downside is that such heavy particles 
Ni cannot be studied in the laboratory. Out of the 7n — 3 parameters in F and Mm, one can 
under ideal conditions probe the 3 neutrino masses and 6 angles and phases in Uj^ experimentally 
(being very optimistic in case of the Majorana phases). The generated asymmetry r]B (the only 
observable number in leptogenesis) can in general depend on all 7n — 3 parameters. Thus, the 
perspectives to constrain A'^z-properties from leptogenesis are very limited. We therefore only 
recapitulate the basic ideas and refer the interested reader to the reviews named above. 

Minimal "Vanilla" leptogenesis - In the simplest scenario one assumes that the Nj are 
very heavy, have hierarchical masses {Tew ^ Af/>i) and the washout is strong. The 

latter statement is usually parametrized in terms of the parameter Kj = T i\t=q/ H\t=Mj — 
{F^F)iiv^/{Mi X lO^^eV), where Tj is the thermal width of Nj particles, H ~ IM^T'^/Mp 
the Hubble rate and g^, the effective number of relativistic degrees of freedom in the plasma 
{g^ = 106.75 in the SM at T > Tew)- If (i^^i^)ii ~ eV x Mi/v"^, as suggested by the seesaw 
relation (jl3p . then Ki ^ 1 and one is in the strong washout regime. In this case the BAU 
responsible for the observed rjs may have been created by the freezeout and decay of the lightest 
sterile neutrino A'^i alone. This is possible if any pre-existing asymmetries, the asymmetries 
generated during the production of the Nj and the asymmetries from the A'^/>i-decays are all 
washed out efficiently by iVi0. This scenario has the advantage that tjb is essentially independent 
of the initial conditions. A major contribution to the final asymmetry comes from the decay 
diagrams shown in figure [TOl In addition, there are lepton number violating scatterings. However, 
vanil la lep togenesis can be qualitatively understood by considering t he d ecays and inverse decays 
only [36j]. This mechanism can reproduce r]B for Mi > 4 x 10^ GeV which requires a rather 



large reheating temperature if Nj are mainly produced thermally ("thermal leptogenesis"). There 
is also a constraint on the mass of the lightest active neutrino as mi < 0.1 eV [366, 363]; see 368 1 



for an review on the connection with neutrino masses. In the weak washout regime Ki < the 
predictive power is much smaller because asymmetries from processes at earlier time (nonthermal 
production during reheating, thermal Nj production, Njyi freezeout and decay...) are not washed 
out efficiently and contribute to tjb- 

Flavour effects - In the simplest scenario it is assumed that only one sterile neutrino A^i 
dynamically participates in leptogenesis and all three active flavours can be treated equivalently 
("unflavoured regime"). In the unflavoured regime, only the linear combination ~ Fai\lLa) of 



''^Though there are no A'^7>i particles in the plasma at T -C Mi>i, the existence of more than one RH neutrino 
is crucial to provide CP- violation in the lepton sector. 

■'^This is possible in spite of (|13|l because Mm and F are matrices. 
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active leptons that couples to A^i is relevant for leptogenesis. In the corresponding flavour basis 
one can ignore the directions in flavour space perpendicular to that, and the problem is equivalent 
to the one flavour case. Both of these assumptions do not hold in general. 

Active flavours have to be treated separately when the Hubble rate drops below the rate at 
which interactions that distinguish active flavours, mediated by the charged Yukawa couplings, 
occur ("flavoured regime"). These interactions destroy the coherence of the flavour state that 
couples to A^^i because they have different strength for the different lL,a- For the r- Yu kawa , this 
happens below T ~ 10^^ GeV. The importance of fl avou r has been realized in [ieg'] and [stJH] 
and meanwhile been studied by various authors, see 362l . l363j | for details and references. Once two 
flavour states are distinguishable, one has to treat the asyi nmet ries stored in each individually. 
They can differ considerable and even have opposite signs 3741 ]. This can affect r]B if e.g. the 
washout is very different for different flavours; it makes leptogenesis possible even if the source 
term does not violate Yla ^a- Flavour e ffects can reduce the lower bound on Mi in generic seesaw 



scenarios by 1 — 2 orders of magnitude 356l . 13751 ] . which is still far out of experimental reach. 

Also the assumption that the generation of today's BAU only involved A'^i-dynamics does 
not hold in general. In the unflavoured regime, A'^i can only wash out the asymmetry that 
is stored in the combination FaiYa in flavour space to which it couples. This singles out a 
particular direction in flavour space. If pre-existing asymmetries or thos e produced by A^7>i have 
a component orthogonal to this, A^i cannot wash them out efficiently 376l . 13771 ]. In fact, it is 



a rather special case that there is no such component. The effect of active flavours in scenarios 
where the heavier Nj contribute has been discussed in |378l4382] . 

Flavour effects offer ways to circumvent the lower bound on Mj. One possibility is that the 
BAU originates from a purely flavoured asymmetry; i.e. the total asymmetry is vanishing or small 
(L ~ 0), but the left handed asymmetry is non- vanishing {\Ll\ ^ |L|). Since sphalerons only 
couple to LH fields, they convert part of into a LH baryon asymmetry. This is possible for 
Mj < V because in this case the seesaw relation ()13p enforces Yukawa couplings that are so small 
that L]^ and Lr m ay n ot equilibrate before sphaleron freezeout. This makes leptogenesis possible 
even for Mm = [383*] (Dirac leptogenesis), and it is also the basis of scenarios discussed in the 
following section 15.21 The mass bound may also be circumvented if th e ge nerated asymmetry is 
resonantly enhanced by a degeneracy between sterile neutrin o rnasses 3841] (re sonant leptogene- 

|385l. lisej: 



sisj. This enhancement is discussed from first principles in [3851, 1386[; in [385|] it is shown that 
the maximal enhancement (compared to the non-resonant case) of the total asymmetry L can be 
expressed in terms of the parameteiEfl 



M1M2 



2(Miri + MaFa) 



(30) 



The enhancement allows values of Mm in the TeV range, which raises hope to probe Ni in high 
energy experiments. In this case it is of course not only the lightest sterile neutrino A^i that 
generates the asymmetry, but the interplay between the two mass-degenerate states, and flavour 
effects in the sterile sector have to be taken into account. This is generally the case when the 
spectrum of Mj is not hierarchical (|Mi — M2I < Mi), not only in the extreme case Mi ~ M2. 
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Figure 11: Constraints on the sterile neutrino masses and mixing = tr(6lt6l) for n = 2 from 
baryogenesis; upper panel - normal hierarchy, lower panel - inverted hierarchy. For the displayed 
case of n = 2 RH neutrinos, baryogenesis can only be successful if their masses are degenerate, the 
para meter M refers to their mean value. For n > 3 RH neutrinos no mass degeneracy is required 
\392ij . The observed BAU can be generated in the region between the solid blue "BAU" lines. The 
regions below the solid black "seesaw" line and dashed black "BBN" line are excluded by neutrino 
oscillation experiments and BBN, respectively. The areas above the green lines of different shade are 
excluded by direct search experiments, as indicated in the plot. The solid lines are exclusion plots for 
all choices of parameters, for the dashed lines t he ph ases were chosen to maximize the asymmetry, 
consistent with the blue lines. Plots taken from [233 /. 
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5.2 Leptogenesis during Nj production 

If some masses M/ are near the electroweak scale or below, the seesaw relation (|13p requires the 
entries F^i to be very small to be consistent with bounds on the active neutrino masses. That 
means that the rate Tj ~ {F'^F)jjT of thermal Nj production is so small that these particles 
may not come into thermal equilibrium until T ~ Tew- During the thermal production before 
equilibration, the nonequilibrium condition is fulfilled and asymmetries La are generated. The 
total asymmetry L is suppressed by Mj/T unless there is a resonant enhancement a la (jSOp . but 



the flavoured asymmetries are unsuppressed [3921 ]. In contrast to the Mj ^> Tew scenarios 
discussed in section 15.11 L is not conserved near Tew j once the Nj reach thermal equilibrium, 
the La are washed out. A baryon asymmetry can be generated despite L ~ if the Nj have not 
thermalized or this washout is not efficient enough to eliminate Ll before sphaleron freezeout 
(because sphalerons only see L 7-, ~ —Lfi 7^ 0). This scenario is often referred to as bary agenesis 



via neutrino oscillations 393l |. though oscillations are not always crucial. One advantage of this 



mechanism is that the Ni are light enough to be within reach of laboratory experiments. This 
makes it one of the few known models of baryogenesis that are empirically testable (the other 
two much studied testable mechanisms are resonant leptogenesis and electroweak baryogenesis 



3941 ]). The perspectives to study the Nr r esponsible for baryogenesis in the laboratory have 



been studied in detail in [23d . I231l . I392l . 13951 ] . If the Nj interact only via the Yukawa couplings 
F, then a clear detection is realistically only possible for masses Mj < a few GeV unless F has 
some spe cial structure leading to cancellations in (jl3p that allow for larger individual elements F^i 



158l . [l59l ]. If they have additional interactions, masses up to a few TeV are in reach of high energy 



experiments, see section \3M and references therein for more details. In the minimal scenario with 
n = 2 sterile neutrinos a mass degeneracy ]Mi — M2\/{Mi + M2) ~ 10~^ is required to produce 
the observed BAU (or 775). Bounds on mass and mixing from the baryogenesis requirement in 
this case are shown in figure [TT] together with other experimental and astr ophysical constraints. 
For n> 2 and Mj in the GeV range no mass degeneracy is necessary 392l |. 



The number of possible leptogenesis scenarios increases greatly when ([T]) is extended by ad- 
ditional degrees of freedom or embedded into a bigger framework (e.g. supe rsyrametr y) ; we do 



not discuss these here in detail and refer the interested reader to the reviews |36ll4363l ]. 



Finally, the CP- and B-violation that make leptogenesis possible can also "work backwards" 
and rule out or constrain other baryogenesis scenarios because they wash out matter-antimatter 
asymmetries created by those mechanisms at higher temperatures. The mass spectrum of the Ni 
can be constra ined by the requirement that they do not wash out the baryon asymmetry in the 
early universe though these conclusions rely on assumptions about additional interactions 



of Nj. A similar argument was suggested in |397 ] to constrain the properties of a fourth neutrino 
generation. 

5.3 Towards a quantitative treatment 

Transport equations - Most quantitative studies of leptogenesis solve a set of momentum 
integrated Boltzmann equations ("rate equations") to predict ijb as a function of the parameters 



''^See [384 l385l . l387H39l| for earlier discussions of this point. 
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in ([T]). In the case of "vanilla leptogenesis" , these can be written as j364i | 

HX^ = -ri(yi-yn (si) 

HX'^^^ = eiTi{Y^-Yl'')-cwTiYB-L. (32) 

Here Yi is the abundance of A^i particles, i.e. the momentum integral over the phase space 
distribution function divided by the entropy density s. Y^'' is its value in thermal equilibrium 
and Y^-L the difference between baryon and lepton abundances (particles minus antiparticles 
for each of them). We use the variable X = Mi/T instead of time, which leads to a factor 
dX/dt ~ HX on the left hand sidj^. ei is a parameter that characterizes the amount of CP- 
violation, 

^ 167r(FtF)ii^ V Jii Ml l&^{F^F)uv^ ^ Aij ^ ) 



The constant cw in the simplest case (washout by inverse decays) is given by 1 /2 times the ratio 
between the number densities of A'^i and active leptons. Calculation of the final Yqt requires 
solving the rate equations ([3T]) and ([32]) . it can be estimated as Yb-l ~ i^^i/ 9* 355 ] £3 where 
the number k < 1 is called efficiency factor and is related to the washout. 

Equations ()3ip and (j32p can be used when only one sterile neutrino dynamically contributes 
to the BAU and the SM leptons can be described in an effective one flavour treatment (un- 
flavoured regime). In the fully flavoured regime flavour dependent interactions act sufficiently 
fast (compared to the time scale ~ l/Pi related to the A'^i dynamics) to fully decohere the dif- 
ferent contributions in the state -FqiI/lo) that couples to A^i. Then (j32]) has to be replaced by 
three different equations for the asymmetries Y^ in the individual flavours. In the intermediate 
regime between these two cases, there is no full decoherence. In this case flavour oscilla tions have 
to be taken into account. This is usually done by density matrix equations |398| . see 369 . 372 1 



for early applications to leptogenesis. These are matrix valued generalizations of the rate equa- 
tions ()3ip and ()32p in which the lepton charges Y^p carry two indices. The diagonal elements 
Yaa of these are simply describe the lepton abundance in flavour a, the off-diagonal components 
describe correlations between different flavours. A flrst principles approach that al lows to study 



the flavoured, unflavoured and intermediate regimes consistently was presented in 399l |. 

If more than one RH neutrino is relevant for leptogenesis, also the different sterile flavours have 
to be treated independently. For a hierarchical mass spectrum, this usually amounts to simply 
replacing (|3ip by n different equations for the abundances Y[ of all relevant sterile flavours Ni. 
In scenarios of resonant leptogenesis or leptogenesis during Ni production, oscillations between 
the different flavours may be relevant. In these scenarios, leptogenesis typically happens at 
temperatures T ^ M/. In this situation transitions between the different helicity states are 
strongly suppressed and the different helicity states of Ni evolve independently. Effectively, they 
act as "particle" and "antiparticle" for the Ni, though this notion can of course not be taken 

521 



literally for a Majorana field I. Correlations between different helicities are usually negligible. 



^°The expression for dX/dt can be complicated when <;« changes during leptogenesis. This can e.g. afltect the 
generation of lepton asymmetries at T <C Tew - 



Our notation is more close to that in [36. 
^At T ^ Mi the average momentum of pf 
hence vr effectively behaves like at massless Weyl field up to corrections 0[Mm/T] 



^■^At T ^ Mi the average momentum of particles is so large that the Majorana mass is kinematically negligible, 
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and it is justified to describe the Nj by two matrices Yn and fo^" the two hehcity states. In the 
mass basis, the abundances Yj can be identified with the elements (liv)//- In the fully flavoured 
regime, where leptogenesis during Nj production and resonant leptogenesis usually take place, 
the kinetic equations then read 



iHX 



iHX- 



dYM 
dX 
dYf, 



iHX 



dX 
dX 



[Hn,Yn\ - -{Tn, Yn - r^"} + -Y^T% , 



I U* V-] " /F* V- V^1\ " V T"^* 



iTlY^ + itv tliYr,-Y;,') 



itr 



N 



(34) 
(35) 
(36) 



The flavour-matrix is the dispersive part of the effective A'^/-Hamiltonian, which leads to 
sterile neutrino oscillations; the rate-matrices F^v, Fjv and form the dissipative part, which 
acts as collision term. These transport coefficients have to be computed from the real- and 
imaginary p arts of the Ni and 1^^^ self energies in thermal field theory. More precise definitions 
are given in 



231[, along with a derivation of ([M|) -([36 |) . 



The equations (j3T]) - (j36]l are rate equations for momentum integrated abundances. They pro- 
vide a good approximation if the momentum distributions are proportional for equilibrium dis- 
tributions (this is often called kinetic equilibrium). If this is not th e case, ea ch momentum mode 



has to be tracked independently by a Boltzmann equation, see e.g. 400|, |40l[. In numerical simu- 
lations, one of course has to sample a finite number of representative momenta. In addition, one 
needs to calculate the collision term as a function of momentum and temperature. This makes 
the treatment technically much more challenging. 



Conceptual issues - A question that has received much attention in recent years is whether 
Boltzmann equations for particle abundances are in principle suitable to describe leptogenesis. For 
instance, calculations based on Boltzmann equations are plagued with a double-counting problem 



402( 1 . The reason is that Boltzmann equations somewhat artificially distinguish the external and 



internal lines of a Feynman diagram; only the external lines are associated with physical particles 
that appear in the phase space distribution functions. If one simply introduces collision terms for 
A'^z-mediated scatterings in addition to those for Nj decays and inverse decays shown in figure fTOt 
one counts the same processes twice: If the intermediate Nj in a scattering is on-shell, then the 
scattering Zq,$ — )• ^Ip is identical to an inverse decay la^ — )• Nj followed by a decay N[ — )• ^Ijs 
and should not be counted independently. This particular problem can be fixed by hand within 
the Boltzmann approach by performing a real intermediate state subtraction (RIS) and using 
finite temperature field theory to calculate the amplitudes, see 403| and references therein for a 
detailed account of this issue. 

In leptogenesis, the leading order processes that contribute to the generation of a matter- 
antimatter asymmetry comes from an interference with loop diagrams as e.g. shown in figure 
[TOl hence leptogenesis is a pure quantum effect. This is in contrast to many other processes 
that are well-described by Boltzmann equations, such as CMB decoupling or BBN. Boltzmann 
equations are semi-classical. The dynamical quantities are phase space distribution functions for 
particles and antiparticles, i.e. classical quantities. The collision terms, on the other hand, are 
calculated from S-matrix elements. This treatment is based on a number of assumptions that 
may be questionable in the dense primordial plasma. The definition of particle numbers and 
the S-matrix are both based on the notion of asymptotic states, the meaning of which is not 
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clear in a dense plasma, where particles are never "far away" from their neighbours. Even if the 
plasma can be described as an ensemble of (quasi)particles, the dispersion relations of these differ 
considerably from those in vacuum. Finally, the collision terms are affected by thermodynamic 
effects (such as Bose enhancement, Pauli blocking. Landau- Pomeranchuk-Migdal effect, possible 
enhancements from multiple scatterings...) and cannot be calculated from the (vacuum) S-matrix. 
The range of validity of the Boltzmann equations and size of possible corrections cannot be 
estimated within this framework and requires a fir st principles tr eatment. This has lead to a great 
interest in nonequi l ibrium quantum fi eld theory 259l . 140444081 ] and applications to leptogenesis 
2061 . [isgl . [ssi, H, [sgi, BH ] in recent years. 



At a fundamental level, the state of a quantum system can be described by an infinite tower 
of n-point correlation functions. In practice, it is usually sufficient to consider two-point func- 
tions. The leptonic charge can conveniently be described in terms of the correlation functiorl^ 
S^{xi,X2)ai3 = \(lL,a{xi)lL,j3{x2) - lL,i3{x2)lL,a{xi)) , where the average (...)= Tx{q . . .) includes 
quantum mechanical and thermodynamic fluctuations (f? is the density operator). The definition 
of this function does not depend on any notion of asymptotic states or (quasi)particles. Prom the 
Fourier transform of this function in the relative coordinate xi — X2 one can define the momentum 
integrated lepton abundance matrices 



Yiit) = 



d^p tr[705+(p;t)] 
(27r)4 s 



dp 



p2 fdpotTb'^Stip;t)] 



27r2 / 2tt 



(37) 



Here Xi = (tj,Xj), t = {ti + t2)/2LJ and the trace runs over Dirac indices. The diagonal elements 
Yiit)aa correspond the lepton numbers in flavour a, the off-diagonals describe correlations between 
different flavours, related to flavour oscillations. For the sterile neutrinos we define the propagator 
G~j^j{xi,X2) = ^{Nj{xi)Nj{x2) — Nj{x2)Nj{xi)). The matrices Yat and can be extracted from 
it as 



YN{t) 



d^p tr[P+7°G+(p;t)] 
(27r)4 s 



d^p tr[P_7°G'+(p;t)] 
(27r)4 ~s 



(38) 



where P± are helicity projectors. The time evolution of and is governed by the Kadanoff- 
Bayi n equation s j424] . which can be obtained from a Dyson-Schwinger equation on a complex time 
path liia - Gii ]. The Kadanoff-Baym equations are exact integro-differential e quations; There are 
different ways to derive effective kinetic equations for Yn, Y^ and Yi from them 259 . 386 . 3991 . 403, 
4171 . l42ll ]. In the fully fiavoured or unflavoured regime, one can neglect active flavour oscillations 
and flnd a flavour basis in which Yi is diagonal. In this case we can, assuming that sphalerons act 
rapidly on the time scale related to the A'^z-evolution with a rate that is in good approximation 
given by the equilibrium onj^. relate Yi to the quantities in ([5^ - ([5S|) as Yq, = (Y5)q,q, — (ns— n^)/3. 

Though the above approach provides a controlled approximation scheme to formally obtain 
effective kinetic equations from first principles, it is not yet a complete theory of le ptogenesis 
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233|, 



because the computation of the transport coefficients in practice can be co niplicated . 

431 4436 ] . It requires knowledge of the quasiparticl e spectru r n in theplasma |2oi, GHl [3- 
439( 1 and inclusion of all processes in the plasma O, H H, Q. In the regime Mj <C T 
it is still not clear which effect corrections from soft and collinear gauge boson exchange have 

'ioTl. liii - tis^ . 



''^We use the notation of [417 
^'*In t he ho mogeneous and isotropic early universe there is no dependence on xi + X2. 
See I429I . I43Q] for some discussion of this point. 
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T he c onclusion on this point is that it is p ossib l e to describe leptogenesis in the unflavoured 

lil?! ]. flavoured |39^, H and resonant [ssi, [sS^] regimes by effective kin etic equations. 



By that, we mean diiTerential equations that are local and of first order in time |259l |. which is 
in contrast to the non-local second order fundamental Kadanoff-Baym equations. The situation 
is relatively simple if the lepton asymmetry is generated at T <C Mj and t he Mj sp ectrum is 
hierarchic al. T hen first principles calculations suggest that in the unflavoured 4141 . |417I | and fully 
flavoured 399f | cases standard Boltzmann equations give results which are correct up to factors 
0[1], provided that the RIS is performed consistently. When flavour oscillations are relevant, one 
has t o use the ir matrix-valued generalisation, the density matrix equations, as e.g. derived in 
|23ll . H l42ll ]. For T > M or when two masses are degenerate, one can still formulate effective 
kinetic equations, such as (I34p -()36p. However, inserting S-matrix elements into the classical 
Boltzmann equation may not give a correct quantitative description. Instead, they should be 
derived using controlled approximations to the full nonequilibrium field theory, which take into 
account the modified (quasi)particle spectrum and dispersion relations in the plasma, the effect 
of (possibly multiple) scatterings with quanta from the therm al bath a i id o t her t hermody i iami c 
effects. T hough great progress has been made in this regime 2071 . l23ll . l385l . l386l . |392| . |415| . 1434 - 
43^, H, the accurate quantitative description is not complete at this stage. 



6 RH Neutrinos as dark matter 

Over the past 80 vearj^. overwhelming evidence has accumulated that most of the mass in 
the observable universe is not composed of baryonic matte^l as we know it. The presence of 
large amounts of non-luminous matter can on one hand be inferred by comparing the observed 
gravitational potential^ to the density of visible matter. Independently of that, it is also required 
in order to explain the clustering of matter and formation of structures in the universe, e.g. 
the power spectrum shown in figure El This process can be studied by the analysis of density 
perturbations in the CMeI^ and by observations of the structure in the universe, such as galaxy 
surveys, gravitational lensing or the absorption in the spectra of distant quasars ("Lya forest"). 

Numerous attempts to simultaneously explain these phenomena by modifications of gravity 
or the presence of compact macroscopic objects (such as lonely planets, black holes or other 
non-luminous star remnants) have failed. In contrast to that, they can easily be understood if 
one assumes the existence of (one or several) new particles that are massive, electrically neutral, 
long lived (compared to the age of the universe), collisionless and have a free streaming length in 
the matter dominated era that is sufficiently small to be consistent with the observed structure 
in the universe. 

One can qualitatively distinguish three types of DM candidates. Cold Dark Matter (CDM) is 
composed of particles that were nonrelativistic at the time of decoupling. Hot Dark Matter (HDM) 

^^The term "dark matter" was already used in [44 ll . 144^ ]. but in that context simply referred to non-luminous 
matter. Zwicky concluded that the velocities in the Coma-system are an "unsolved problem", which, however, did 
not receive much attention for decades. 

^^In this context all matter that is composed of SM fermions (including leptons) is usually referred to as "bary- 
onic" . 

^*The gravitational potential can be studied by tracking the dynamics of astrophysical objects on various scales, 
such as stars in galaxies (in particular rotation curves) as well as galaxies and ionized gas clouds within galaxy 
clusters, and by gravitational lensing (weak and strong). 

^^The fact that most of the matter started to cluster before the decoupling of photons implies that it had 
decoupled considerable earlier, hence is not composed of (electrically charged) SM particles. 
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particles were relativistic at the time of decoupling and remain so into the matter dominated 
epoch, when structures can grow nonlinearly due to gravitational collapse. Warm Dark Matter is 
relativistic at the time of freezeout, but becomes nonrelativistic during the radiation dominated 
epoch. 

CDM scenarios p redict tha t smaller scale objects form first and then merge into bigger struc- 
tures ("bottom up") [443 . 444 1 . This is in good agreement with observations as we ll as r ium erical 
simulations of LSS. On smaller scales, it was noticed more than a decade ago 4451 . 44d | and 
remains an unsolved puzzle that CDM simulations tend to predict more objects (satellite galax- 
ies and subhalos) than observed. Furthermore, they do not reproduce shape of galactic halos 
("cusp/core issue"). However, at this stage this is far from being a reason to disregard CDM 
because the discrepancy may be related to the resolution or other technical limitations of the 
simulations (to date, most of them are pure DM simulations that do not include baryons consis- 
tently). Furtherm ore, small objects are hard to observe if they fail to confine gas and form stars, 



see e.g. 



447l4450i . thou gh th is may not explain the discrepancy for objects that are too big to fail 
to attract gas, see e.g. [4511] . A recent summary on issues and possible solutions in the ACDM 
model can be found in j452l |. 

In contrast to that, HDM predicts that large structures form first because primordial fiuctu- 
ations in the gravitational potential have been erased due to the large free streaming length of 
DM particles ("top down"). This is in contradiction to the observed LSS, which rules out the 
only potential DM candidate within the SM, the neutrinos i^lS In WDM scenarios, structure 
formation on scales larger than the particles' average free streaming length Xdm is similar to 
CDM, but differences on smaller scales are expected to be visible in various observables, such as 
the matter power spectrum, halo density profile, halo mass function, subhalo density profile and 
subhalo mass function. Qualitatively, structures on scales smaller than Xdm are suppressed due 
to the free streaming. 

Sterile neutrinos are collisionless and can be very long lived, hence they are an obvious DM 
candidate. This s c enario has been st udied by a large number of authors, see e.g. 



Ill, UM, 



isi . [293 . [333I . [345I . lisi - tiTi . I47i - t493 ] . Formally, sterile neutrinos are WDM candidates, though 



this classification is slightly ambiguous because momentum distribution may be very different 
from Fermi-Dirac. Their properties are constrained by the requirements outlined below. In what 
follows, we require that all DM is composed of sterile neutrinos. If they make up only a fraction 
of r^DM; then the bounds weaken considerably. For instance, sterile neutrinos with eV masses 
could give a subdominant HDM contribution to Qdm if there are other CDM particles that ensure 
consistency with structure formation. 



Stability - Nj particles are unstable and decay via the ^-suppressed weak interaction. If they 
are DM, their lifetime must be longer than the age of the universe. 



X-ray bounds - The radiative decay Nj — )• jy^ via processes as shown in figure [T3l predicts the 
emission of photons with energy Mj /2 from DM dense regions 
of such signal in the data of different X-ray observatories (such as XMM, Chandra and Suzaku) 

^''The active neutrino background can be viewed as a small HDM contribution to the matter density in the 
universe. 

^^Here we assume that i/_r have (at least) the Yukawa interactions in ([T}. If such interactions are suppressed, the 
RH neutrinos that constitute DM can be much heavier and may even be related to the observed excess [g] of 7-ray 
emission near the galactic center [49 5^ ). 
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The non-observation 
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Figure 12: Different constrain ts on sterile neutrino DM mass and mixing from X-ray observations. 
The constraints are take from \46{\l (red), rescal ed by a factor of two due to mass estimate uncertain- 
ties as recommended in fJ87l /. \48'A] (grey) and [45?t / (pur ple). Som e analyses have cla imed stronger 
constraints, but were later found to be too optimistic. In [45^. 46j|j it was found that \45d] underes- 
timated the flux by two orders of magnitude. According to \46i ] the mass was overestimated in f46C( / 
leading to too restricti ve c onstraints. The constraints in [48a] might be too restrictive due to the 
choice of source profile \46i] . The spectral resolution seems to be overestimated in NgSl/. cf. Chandra 
Proposers Guide; this seems to be the main reason for the stronger bounds used in 1120]. Thanks to 
Signe Riemer-S0rensen for the plot and comments. 




Figure 13: Example for a contribution to the radiative decay of a DM sterile neutrino Ni 
The coupling of Ni to the W-boson is suppressed by the mixing 6a\ as indicated. 



imposes upper bounds m the mass-mixing plane; 



these are displayed in figure [T2l In figure [HI X-rav bounds are combined with other constraints. 
For decaying DM, the signal strength scales only linear with the DM column density along the 
line of sight, hence one can expect a signal from wider range of astrophysical objects than for 
annihilating DM. This makes future searches promising despite the fact that there are various 
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1 2 5 10 20 50 

Ml [keV] 

Figure 14: Different constraints on sterile neutrino DM mass and mixing, assuming vr have no 
interactions in addition to the Yukawa couplings in 0. The blue region i s excluded by X-ray 
observations, the dark grey region M < I keV by the Tremaine-Gunn bound \47i . \474 . \49d l. The 
solid black lines are "production curves" for thermal production. For all points on the upper black 
line the obse rved ^Idm is produced in the absence of lepton asymmetries (for ~ 0, non-resonant 
production) \23A] . Points on the lower solid black line yield the correct Q,dm for \Ya\ = 1.24 • 10~* at 
T — 100 MeV, the dashed line for \Ya \ — 7- 10~^ . For these values the resonant production mechanism 
contributes. The region between these lines is accessible for intermediate values of Y^ . We do not 
display bounds derived from structure formation because they depend on Ya in a complicated way 
and there are considerable uncertainties. Plot taken from \23ul l 



astrophysical s ources in the keV range. The perspectives for future searches have e.g. been 
outlined in Finally, it is worth to emphasize that this scenario is falsifiable: 

if DM is made of RH neutrinos that are produced thermally in the early universe through their 
mixing with active neutrinos, their mass and mixing are constrained in all directions in figure 
and this scenario can be found or falsified on human time scales. 



Phase space analysis - As a fer mionic DM candidate, the mass of RH neutrinos is constrained 
by the Tr ema i ne-G unn bound 49 d ] on the phase space density in the Milky way's dwarf spheroidal 
galaxies 472 . 474]. This yields a lower bound of Mn > 1 keV. In {490 1 it has been shown that 
phase space arg ume nts on the galactic scale indeed favour a keV DM particle. The same authors 
have argued in [49l[ | that, due to the high phase space density, a quantum mechanical treatment 
is necessary on these scales, which may solve the "cusp issue" . 



Production in the early universe It is not known when and how DM was produced in the 
early universe. If it is composed of sterile neutrinos, then there are several possible production 
mechanisms. 

• Thermal production via mixing (non-resonant) - If i^/j have nonzero Yukawa couplings 
F (i.e. rriD ^ 0), then they are produced thermally from the primordial plasma via their 
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mixing with the SM neutrinos vl 493] • If this production was efficient enough to bring 
them into thermal equihbrium, their density would be bigger than VLdm and "overdos e" the 
universe unless they are diluted by entropy production at some later stage jioi - lsoot - This 
puts an upper bound on the mixing angle, given by the upper production curve in figure [111 
The production can yield the correct ri/jjv/ if ^ is small enough that the DM sterile neutrino 
never reach thermal equilibrium. In that c ase they still have a momentum distribution that 
is proportional to a Fermi-Dirac spectrum 454 . 497 . 501 . 502| . which makes them a WDM 
candidate that is at least disfavoured by structure formation arguments (see below). 

Resonant thermal production - The properties of (q uasi)parti cles in the primordial 
plasma are modified by the interactions with the medium [SOSySCyl. In the presence of a 

' I ca n lead to a level 

crossing between active and sterile neutrino dispersion relations |233l . l508l |. This results in 
a resonant enhancement of the sterile neutrino production rate. The res ultin g spectrum is 
nonthermal, with higher occupation numbers for low momentum modes 2331 ] . Production 
curves for two different values of the asymmetry are also given in figure [TH Note that these 



have to be taken with some care because in the calculation of the production efficiency [2331 ] 



it was assumed that Y^, = Yjj = Yt-. In reality the a symm etries in different flavours may be 
different and can even have opposite sign (see e.g. [23]]] for a particular model). 

Thermal production beyond the SM - The fields vn are singlet under SM the gauge 
group and only interact with SM particle s via their Yukawa couplings F. However, if 
there exists and extended "hidden sector" 509| or extended Higgs sector, they may have 
additional interactions with these. They may also be charged under some extended gaug e 
group that is broken at high energies, as e.g. in the left-right-symme tric model 482 ]. 
Such additional interactions would contribute to the thermal production [499^ in the early 
universe. Since they usually increase the production rate, the universe may ove rdos e unless 
the DM abundance is diluted by additional entropy production at later times [5101 ]. 



Nonthermal production - RH neutrinos can al so be produced non therrnally , 
a coupling to an inflaton 511145131 ] . the SM Higgs 514], other scalars 15lU515l | 
gravity (si 



e.g. due to 
or modifled 



All these scenarios are constrained by the requirements to produce the correct ^dm and a mo- 
mentum distribution that is consistent with the observed LSS. 



Structure formation - The masses dictated by the above bounds suggest that, if RH neutrinos 
mix with LH neutrinos, then they formally are WDM, but their momentum distribution can be 
non-thermal. Structure formation in WDM scenarios on scales above the free streaming length 
is similar to CDM. On smaller scales the formation of structures is affected in different ways, the 
study of which is in general complicated due to the nonlinear nature of the clustering process. 
Most importantly , ori e expects a cut in the matter power spectrum at scales below the free 
streaming length 5171 ]. In addition, the suppression of small scale structures should be visible 
in the halo mass func tion (which counts the number of haloes per unit volume per unit mass 
at giv en redshif t) 518| and affect the gravitational collapse leading to the formation of the first 
stars [520I ]. Unfortunately, these arguments do not directly constrain M^r, but the free 
streaming length. Furthermore, structure formation is a highly nonlinear problem in the regime 
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5pk/ p > 1 that can o nly be studied quantitatively with expensive many body simulations. With 
few exceptions 292], the se s imulations assume that the DM distribution is proportional to a 
Fermi-Dirac distribution 454l |. I f this is th e case, then the comoving free streaming length Xdm 
is simply related to the mass by [494l . l517l | 



A 



DM 



IMpc 



keV\ (pn) 



M 



N 



(39) 



where pN and py are the spatial momenta of sterile and active neutrinos. Such spectra are 
produced by the non-resonant thermal production from the pr imordial plasma. In that cas e Lyg 
forest observations strongly constrain the viability of WDM I46i4466l . lirol . l479l . lisol . [521 



52d |: for sterile neutrinos they impose a lower bound of > 8 keV 47Clll . The n it may be very 



difficult to distinguish RH neutrino WDM from CDM observationally 470|, l527l | . If, on the other 
hand, RH neutrino DM is produced by some other mechanism and has a nonthermal spectrum, 
the relation (j39p between mass and free streaming length does not apply, and simulations of 
structure formation based on thermal WDM do not allow to draw any general conclusions. In 
case of the resonai it therma l production mechanism, the resulting spectrum has been calculated 
for Ye = Yf, = Yr lis^, H . It can be viewed as a superposition of a WDM component and a 
nonthermal cold component [470|,|47l|, where the Lya forest bounds allow thewarm component to 
be much "warmer" than in conventional WDM scenarios. First simulations 292l | indicate that this 
scenario may perform better than CDM in predicting the abundance of small scale structures (such 
as satellite/dwarf galaxies), but the question is not settled at this stage. On one hand, only little 
is known about how much structure actually exists at sm all (s ub-galactic) scales. While lensing 
flux and stellar streams suggest the existence of subhalos 452], there is no direct obser vation of 
such structures and their existence is disfavoured by some stability considerations 528l |. On the 
other hand there are no systematic studies of structure formation in the nonlinear regime that 
include the nonthermal component. Arguments against sterile neutrino DM are essentially based 
on extrapolations of results from simulations that assume purely thermal spectra. 



7 A theory of almost everything 

Right handed neutrinos vr, described by the Lagrangian ([1]), provide plausible explanations for 
the phenomena (I), (II), (III), (i) and (ii) named in the introduction. In this section we address 
the question how many of these phenomena can be explained simultaneously by RH neutrinos 
alone. 

Minimal case n = 3 - T he Lagrangian ([1]) with n = 3 RH neutrinos and Mm at or below the 
electroweak scale |529l . 53Cll ] is an extension of the SM motivated by the principle of minimality 



(or Ockham's razor); it makes only minimal modifications to the known principles and particles 
in nature. There is no modification to the SM gauge group, the number of fermion families is 
unchanged, there is no modification to the bosonic field content of the SM or the mechanism of 
elect roweak symmetry breaking, and no new scale above the electroweak scale is introduced^. In 
it has been shown that this model can indeed simultaneously explain (I)-(III); that is, all 



known evidence for particle physics beyond the SM. This minimal scenario is known as Neutrino 
Minimal Standard Model (z^MSM). Various aspects of this model have been studied in the past 



This implies that the hierarchy problem of the SM can be avoided [53lH53i 
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by different author s Eli, flij. B E B Esi-flS? .B liii . liij . IsQ^ . H, Sill, [il^, Ssl, 
49^ . [Hoi, [Hnl, I5li.l514l. l529l . l53nl . 15321454^ 866 [l87l.l23lt for a detailed introduction and review. 



In the i^MSM, two sterile neutrinos (^''2,3) have degenerate masses between roughly a GeV 
and the electroweak scale; the third one (A^i) has a mass in the keV range and acts as DM 
candidate. generate active neutrino masses via the seesaw mechanism; at the same time, 

the CP-violating oscillations during their thermal production at temperatures T > Tew produce 
the baryon asymmetry of the universe via the mechanism outlined in section 15.21 The lepton 
asymmetries generated during iV2,3 production are washed out after sphaleron freezeout at T ~ 
Tew- The freezeout (T ~ M2,3 ) and decay (T < M2,3) of the A'^2,3 particles produce new, late 
time lepton asymmetries, which can be orders of magnitude bigger than the baryon asymmetry. 
The late time asymmetries are capable of enhancing the rate of A'^i production sufficiently to 
explain the observed JIda/ for ^ 10~^, while being in agreement with all known constraints 
on DM properties. 

The minimal extension of the particle content in the i^MSM comes at the price of a "fine 
tuning" in the mass splitting of two RH neutrinos: M2 and M3 have to be equ al to a le vel of 
precision ~ 10~^^ (23lli | . This degeneracy may be related to new symmetries 1581 . 5401 ]. but 
cannot be explained in the framework of the i/MSM itself. The tuning mainly arises from the 
requirement to produce a late time lepton asymmetry > 10~^; this is necessary to sufficiently 
enhance the A'^i production to explain the observed i^DM- If so me of the asymmetries generated 
before sphaleron freezeout can be preserved from washout 547l | or A'^i is produced by some other 
mechanism (e.g. during reheating), then the required degree of degeneracy reduces to ~ 10~^. 
This weaker degeneracy is necessary for baryogenesis, i.e. to explain the observed fi^ from CP- 
violating oscillations of two sterile neutrinos A'^2,3- The third sterile neutrino A'^i cannot contribute 
significantly to leptogenesis if it is DM because in this case its couplings F^i are constrained to 
be tiny by X-ray bounds on its decay width, cf. section [6l For the same reason it also cannot 
make a significant contribution to neutrino masses via the seesaw mechanism (I13p : this implies 
that the lightest active neutrino is eff'ectively massless in the z/MSM. 

The particles N2 and can be searched for in collider experiments, cf. section 13.41 Known 
bounds on their mass and mixing are summarized in figure [T5j The DM candidate A'^i is too 
weakly coupled to be studied at colliders, but can be found by indirect DM searches for the X-ray 
emission line from the decay A^i — )• 7V, cf. figure [T3l Its properties are also constrained by phase 
space and structure formation considerations. Figure [T4l summarizes different bounds on A'^i mass 
and mixing. 



More than three sterile neutrinos - The required amount of "fine tuning" can be con- 
siderably reduced if there are more than three sterile neutrinos. For instance, if three sterile 
neutrinos participate in baryogenesis, then there is no need for degenerate masses even if they 
are as light as a few GeV [392]. If one at the same time requires one sterile neutrino with keV 
mass to compose DM, then that means there must be at least n = 4 of them altogether. If one, 
in addition to the evidence (I)-(III), also wants to address the hints (i) and (ii) in this framework, 
then there must be least one more sterile neutrino with an eV mass. This light sterile neutrino, 
responsible for the oscillation anomalies, would leave the thermal history of the universe before 
BBN essentially unaffected due to the smallness of its couplings enforced by the seesaw relation 
P3p . At the same time it can be responsible for the iVeff > 3 inferred from BBN and the CMB, 
see section [H The late time lepton asymmetries generated by the freezeout and decay of its 
Mj > 1 GeV siblings could help to suppress the efficiency of its thermal production and ease the 
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Figure 15; Constraints on the N2,3 masses M2,3 — M and mixing — tr{6^6) in the uMSM; upper 
panel - normal hierarchy, lower panel - inverted hierarchy. In the region between the solid blue 
"BAU" lines, the observed fis can be generated. The observed Qdm can be explained in terms of 
sterile neutrinos inside the solid red "DM" line (there the lepton asymmetry at T — 100 MeV can 
be large enough to resonantly enhance the Ni production). The CP-violating phases were chosen 
to maximize the asymmetry at T = 100 MeV. The regions below the solid black "seesaw" line and 
dashed black "BBN" line are excluded by neutrino oscillation experiments and BBN, respectively. 
The areas above the green lines of different shade are excluded by direct search experiments, as 
indicated in the plot. The solid lines are exclusion plots for all choices of uMSM parameters, for the 
dashed lines the phases were chosen to maximize the late time asymmetry, consistent with th e red 
line. The contributions from Nj to QvPjB-decays in the vM SM are negligible ',116,] or negative \l2!ij 
and do not lead to significant constraints. Plot taken from \23d] . 
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tension between oscillation anomaly and dark radiation signals, see section [^^21 Thus, n > 4 RH 
neutrinos described by the Lagrangian ([T]) alone can be capable of explaining the points (I)- (III) 
as well as (i) and (ii) simultaneously without "new physics" above the electroweak scale; that is 
all confirmed fundamental physics phenomena that cannot be explained in the SM except those 
that are likely to involve accelerated cosmic expansion. Hence, this framework at the current 
stage provides a phenomenologically complete description of particle physics. 

8 Conclusions 

Right handed neutrinos provide plausible explanations for neutrino oscillations, the observed 
dark matter and the baryon asymmetry of the universe, which cannot be understood within the 
SM. They may also be responsible for the anomalies reported by some short baseline neutrino 
experiments and hints for dark radiation in cosmological data. Right handed neutrinos can be 
embedded naturally into bigger frameworks, such as left-right symmetric theories, grand unified 
theories, super symmetric extensions of the SM and various models of neutrino mass generation. 
The different phenomena point towards different right handed neutrino mass ranges, which may 
of course coexist in nature. Several of these can be probed by direct and indirect laboratory 
searches, astrophysical observations and cosmological data in the near future. We summarize 
well-motivated mass ranges and their phenomenological relevance in the overview table appendix 

El 
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A Overview Table: Majorana mass scales and observables 



Mm 


Motivation 


v-oscillations 


laboratory searches 


CMB 


BBN 


DM 


Leptogenesis 


$cV 


"v-oscillations anomalies, 
dark radiation 


massses by seesaw," 
explain anomalies'" 


oscillation anomalies, 
/3-decays 


explain N^g > 


may explain A^eff > S*" 


no 


no 


kcV 


DM 


no if DM<= 


direct searches?'', 
nuclear decays?'' 


act as DM, 
no effect on N^g 


effect on N^g 
too small if DM 


good candidate 


no 


MeV 


testability, why not? 


masses by seesaw 


intensity frontier. 


unaffected 


constrains 
Mj > 200 MeV 


no' 


possible 
(fine tuning) 


GeV 


testability, 
minimality 


masses by seesaw 


intensity frontier, 
EW precision data, 


unaffected 


unaffected 


no' 


possible 


TcV 


minimality, 
testability 


masses by seesaw 


LHC 


unaffected 


unaffected 


no' 


possible 


> lO^GeV 


grand unification, 
"naturally" small v-masses 


masses by seesaw 


too heavy to be found 


unaffected 


unaffected 


no' 


works naturally 



Colour code: green = can affect, red = does not affect 

This assumes that the observed Arrtatm and Amgoi are generated by un other than those responsible for the anomalies, i.e. n > 2. 
' Sterile neutrinos with masses and mixings suggested by the oscillation anomalies would be in thermal equilibrium in the early universe, hence increase A^eff one unit 
per species. In ACDM cosmology, this is in conflict with recent CMB data and IfDM constraints. If some mechanism prevents the sterile neutrinos from thermalizing or 
there are deviations from the standard ACDM model, this conflict can be avoided. 

uji with keV-masses can generate neutrino masses via the seesaw mechanism, but then they are too short lived to be DM. If they are DM, then their Yukawa coupling 
is too small to contribute to the seesaw - they can either be DM or contribute to the seesaw mechanism. 

'' It is disputed whether the signal can be distinguished from the active neutrino background; for the case that keV sterile neutrinos compose all DM the author considers 
it very unlikely that such searches are successful because of the astrophysical constraints on the mixing angle. 

This applies to sterile neutrinos thermally produced via their mixing. Sterile neutrinos with Mj 2> kcV can be DM if _F ~ ensures their stability and the production 
in the early universe is due to an unknown interaction. 



B Dirac and Majorana masses 



This appendix summarizes some basic aspects of Majorana fermions, see [231 . l24l ] for details. 
In the SM, matter is composed of fermions. In relativistic quantum field theory these can be 
described by two component (Weyl) spinors ipL and ipR, which transform under the (irr educ ible) 



left or right handed representations of the Poncaire group, respectively. For historical 5481 ] and 



computational reasons fermions are, however, often represented by four component spinors. In 
this review we adopt this common "overnotation" ; using the Weyl representation of matrices, 
we define four component spinors ^'l = (V'LjO)^ and = (OjV'r)"^- Consider sets of left and 
right handed spinors ^L,i and ^R,j, where the indices j and j run from 1 to n and m, respectively, 
and label individual fields ("flavours"). The most general free Lagrangian reads 

j I 

- (^'L^^L + ^R^^ij) - ^Lmn'i'B - 2 (^L^Af ^'i + "^rMm^r) + h.c, (40) 

where we have suppressed flavour indices {^r, are now flavour vectors with components ^R^i 
and ^L,j)- The matrices niD, rnu and Mm can be combined into a (n + m) x (n + m) matrix 
9Jl, the physical mass squares are the eigenvalues of OJT^OJt. The mass term reads 

i(wr ( ) . 4(*r 15) ( :| ™f ) ( *^ ) . (41) 

All fermions in the SM are charged under some gauge group, thus any intrinsic mass terms are 
forbidden by gauge invariance. Terms as m/j are generated by the Higgs mechanism. For Mm = 0, 
ruM = one can express ()40p in terms of Dirac spinors ^ = ^ r + ^ l as £ = ^{i^ — itid)^- This 
is the reason why one can describe charged leptons and quarks by Dirac spinors though the basic 
building blocks of matter are Weyl fermions: the conserved charge related to the unbroken U{1) 
subgroup in the electroweak theory leads to mass matrix that allows to combine two Weyl fields 
into a Dirac spinor. If, on the other hand, itim^Mm 7^ 0, then one can form n + m Majorana 
spinors [549] . In the simplest case, when m/j = and tum and Mm are diagonal, one can define 
Xi = ^L,i + and ^\>j = ^ R^j + ^^Rj- Formally x and ijj are four component objects, and the 
Lagrangian can be written in analogy to the Dirac equation, i(x(i^— m'M)x+4'(^^— -^m)4')) but 
they obey the additional Majorana conditions x = x'^i ^\) = and have only two independent 
components. Due to this property the neutral fermions can annihilate with themselves (are their 
own antiparticles) . 

Because of these considerations is usually referred to as Dirac mass term, while tum and 
Mm are called Majorana mass terms. However, apart from the two special cases discussed above, 
they have not much to do with the appearance of Dirac or Majorana particles. Even for tud = 
the particle spectrum may contain Dirac spinors. Consider, for instance, a set of three fields ^l, 
^R^i and ^_R,2 with mD^mM = and Mm = Ml2x2- Then one can combine ^_r,i and ^'i?,2 into 
a Dirac spinor W = {i^\>l + 'i[>2)/\/2- On the other hand, one can split any Dirac spinor into two 
Majorana spinors in this way. 

Thus, fermions can generically described by two component Weyl spinors. If one, for compu- 
tational convenience, prefers to use four component spinors, one can use the chiral spinors ^r^l- 
Only in special cases one can combine these into Dirac spinors. This is possible whenever the full 
mass matrix Tl has degenerate eigenvalues (with opposite sign), and not necessarily related to 
the presence of "Dirac" or "Majorana mass terms". 
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